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Abstract 

During batch growth on mixtures of two growth-hmiting substrates, microbes con- 
sume the substrates either sequentially or simultaneously. These growth patterns 
are manifested in all types of bacteria and yeasts. The ubiquity of these growth 
patterns suggests that they are driven by a universal mechanism common to all 
microbial species. In previous work, we showed that a minimal model accounting 
only for enzyme induction and dilution explains the phenotypes observed in batch 
cultures of various wild- type and mutant /recombinant cells. Here, we examine the 
extension of the minimal model to continuous cultures. We show that: (1) Sev- 
eral enzymatic trends, usually attributed to specific regulatory mechanisms such as 
catabolite repression, are completely accounted for by dilution. (2) The bifurcation 
diagram of the minimal model for continuous cultures, which classifies the substrate 
consumption pattern at any given dilution rate and feed concentrations, provides a 
a precise explanation for the empirically observed correlation between the growth 
patterns in batch and continuous cultures. (3) Numerical simulations of the model 
are in excellent agreement with the data. The model captures the variation of the 
steady state substrate concentrations, cell densities, and enzyme levels during the 
single- and mixed-substrate growth of bacteria and yeasts at various dilution rates 
and feed concentrations. (4) This variation is well approximated by simple analyt- 
ical expressions that furnish physical insights into the steady states of continuous 
cultures. Since the minimal model describes the behavior of the cells in the absence 
of any regulatory mechanisms, it provides a framework for rigorously quantitating 
the effect of these mechanisms. We illustrate this by analyzing several data sets from 
the literature. 
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1 Introduction 



The mechanisms of gene regulation are of fundamental importance in biology. 
They play a crucial role in development by determining the fate of isogenic 
embroyonic cells, and there is growing belief that the diversit y of biologica l 



organisms reflects the var i ation of their regulatory mechanisms (ICarroll et al 



20051 : IPtashne and Gand . l2002f ). 



Many of the key principles of gene regulation, such as positive, negative, and 
allosteric control, were discovered by studying the growth of bacteria and 
yeasts in batch cultures containing one or more growth-limiting substrates. 
A model system that played a particularly important role is th e growth of 



Esch erichia coli on lactose and a mixture of lactose and glucose (IMuller-Hill . 



19961). 



It turns out that the enzymes catalyzing the transport and peripheral catabolism 
of lactose, such as lactose permease and /9-galactosidase, are synthesized or 
induced only if lactose is present in the environment. Th e molecular mecha- 



nism of induction was discovered by Monod and coworkers ([Jacob and Monodl . 



19611 ). They showed that the genes encoding the peripheral enzymes for lac- 
tose are contiguous and transcribed sequentially, an arrangement referred to 
as the lac operon. In the absence of lactose, the lac operon is not transcribed 
because the lac repressor is bound to a specific site on the lac operon called 
the operator. This prevents RNA polymerase from attaching to the operon and 
initiating transcription. In the presence of lactose, transcription of lac is trig- 
gered because allolactose, a product of /5-galactosidase, binds to the repressor, 
and renders it incapable of binding to the operator. 

When Escherichia coli is grown on a mixture of glucose and lactose, lac tran- 
scription, and hence, the consumption of lactose, is somehow suppressed until 
glucose is exhausted. During this period of preferential growth on glucose, the 
peripheral enzymes for lactose are diluted to very small levels. The consump- 
tion of lactose begins upon exhaustion of glucose, but only after a certain 
lag during which the lactose enzymes are built up to sufficiently high lev- 
els. Thus, the cells exhibit two exponential growth phases separated by an 
interme diate lag. Mon od referred to this phenomenon as diauxic or "double" 



growth (iMonodl . Il942h 



Two major molecular mechanisms have been proposed to explain the repres- 
sion of lac transcription in the presence of glucose: 



(1) Inducer exclusion (jPostma et al.l . 119931 ): In the presence of glucose, en- 
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zyme IIA^''^, a peripheral enzyme for glucose, is dephosphorylated. The 
dephosphorylated IIA^^'^ inhibits lactose uptake by binding to lactose per- 
mease. This reduces the intracellular concentration of allolactose, and 
hence, the transcription rate of the lac operon. 

Genetic evidence suggests that phosphorylated IIA^'^ activates adenylate 
cyclase, which catalyzes the synthesis of cyclic AMP (cAMP). Since glu- 
cose uptake results in dephosphorylation of IIA^'^, one expects the cAMP 
level to decrease in the presence of glucose. This forms the basis of the 
second mechanism of lac repression. 



(2) cAMP activation (jPtashne and Gannl . 120021 ): It has been observed that 
the recruitment of RNA polymerase to promoter is inefficient unless a 
protein called catabolite activator protein (CAP) is bound to a specific 
CAP site on the promoter. Furthermore, CAP has a low affinity for the 
CAP site, but when bound to cAMP, its affinity increases dramatically. 
The inhibition of lac transcription by glucose is then explained as follows. 
In the presence of lactose alone (i.e., no glucose), the cAMP level is 
high. Hence, CAP becomes cAMP-bound, attaches to the CAP site, and 
promotes transcription by recruiting RNA polymerase. When glucose is 
added to the culture, the cAMP level decreases by the mechanism de- 
scribed above. Consequently, CAP, being cAMP-free, fails to bind to the 
CAP site, and lac transcription is abolished. 



Mechanistic mathe matical models of the glucose-lactose diauxie appeal to both 



19751 : IWong et al. 



1997h 



these mechanisms ( Kreinling et al.l . l200ll : ISantillan et al.l . 120071 : iDedem and Moo- Young 



Over the last few decades, microbial physiologists have accumulated a vast 
amount of data showing that diaux i e growth is ubiquitous (reviewed in lEglil . 
19951 : [Harder and Dijkhuizenl . Il982l : iKovarova-Kovar and Egli Il998l ). It oc- 
curs in diverse microbial species and numerous pairs of suhstitutable sub- 
strates (satisfying identical nutrient requirements). In many of these systems, 
the above mechanisms play no role. For in stance, cAMP is not detectable in 



gram-positive bacteria (IMach et al.l . Il984l ). and has little effect on the tran- 



scriptio n of peripheral enzymes i n pseudomonads (ICoUier et al.l . 119961 ) and 
yeasts ( lEraso and Gancedd . Il984l ). It has been proposed that in these cases, 
the mechanisms of repression are analogous to those of the lac operon, but 
inv olve different tran scription factors, such as CcpA in gram-pos itive bacte- 
ria (IChauyauxl. 119961). Crc in pseudomonads (iMorales et al.l . l2004l ). and MigA 
in yeasts ( I Johnston! . Il999l ). 



Paradoxically, there is growing evidence that cAMP activation and inducer 
exclusion are not sufficient f or explaining lac repression. Not long after the dis- 
cover y of cAMP activation (jPerlman and Pastanl . Il968l : lUllmann and Monodl . 
19681 ) , Ullmann presented several lines of evidence showing that thi s mech- 
anism played a relatively modest role in diauxie growth (jUUmanrJ . Il974l ). 
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But the most compelling evidence was obtained by Aiba and coworkers, who 
showed that the cAMP levels are essentially the same during the first and 
second growth phases; furthermore, glucose-mediated lac repression persists 
in the presence of high (5mM) exogenous cAMP levels, and in mutants of 
E. coli in which the ability o f cAMP to influence lac tr anscription is com- 



pletely abolished (llnada et all Il996l : iKimata et all 119971 ). In Section 13.11 we 
shall show that the positive corre l ation between lac expression and intracel- 



lular cAMP levels (lEpstein et al.l . 119751 ) . which forms the foundation of the 



cAMP activation mechanism, does not reflect a causal relationship: The ob- 
served variation of lac expression is almost entirely due to dilution (rather 
than cAMP activation). 

The persistence of lac repression in cAMP-independent cells ha s led to the hy- 



poth esis that inducer exclusion is the sole cause of repression dKimata et al 



19971 ) . However, this mechanism exerts a relatively mild ef fect. The lactose up- 
take rate is inhibited only ~50% in the presence of glucose (IMcGinnis and Paigenl . 
1969h . 



Thus, transcriptional repression, by itself, cannot explain the several hundred- 
fold repression of lac during the first exponential growth phase of diauxic 
growth. In earlier work, we have shown that complete repress ion is predicted 
by a minimal model accounting for only induction and growth (|Narangl . Il998al . 

2nn6h . 



Yet another phenomenon that warrants an explanation is the empirically ob- 
served correlation between the substrate consumption pattern and the spe- 
cific growth on the individual substrates. For instance, in the case of diauxic 
growth, it has been observed that: 

In most cases, although not invariably, the presence of a substrate per- 
mitting a higher growth ra te prevents the utilization of a second, 'poorer' 
substrate in batch culture ( iHarder and Dijkhuizenl . Il982l ). 



It turns out that sequential consumption of substrates is not the sole growth 
pattern in batch cultures. Monod observed that in several cas es of mixed- 
substrate growth, there was no diauxic lag (IMonodl . Il942l . 119471 ). and subse- 
quent st udies have shown that both substrates are often consumed simulta- 
neously ( lEgli I1995I ). The occurrence of simultaneous substrate consumption 
also appears to correlate with the specific growth rates on the individual sub- 
strates. Based on a comprehensive review of the literature, Egli notes that: 

Especially combinations of substrates that support m edium or l ow maxi- 
mum specific growth rates are utilized simultaneously (lEglil . 119951 ) . 



Recently, we have shown that the minimal model provides a natural explana- 
tion for the foregoing correlations, which hinges upon the fact that the enzyme 
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Figure 1. If the induction rate of the peripheral enzyme is independent of the 
inducer level, it s activity decrease s with the dilution rate, ( a) The activi t ies of 
/3-galactosidase ( Silver and Mateles . 1969I . Fig. 1) and amidase ( Clarke et al. . 1968 . 
Fig. 2) during growth of the constitutive mutants, E. coli B6b2 and P. aerugi- 
nosa Cll, on lactose and acetamide, respectively, (b) The activity of formaldehyde 
dehyd rogenase (FDH) during glucose-limited growth of H. polymorpha and C. boi- 



dinii (jEgli et al 
and (fTSjl . respectively. 



1980 . Fig. lc,d). The curves in (a) and (b) shows the fits to eqs. (fT7 



dilutio n rate is proportional to the specific growth rate (iNarang and Pilyuginl . 
2007al ). Roughly, substrates that support high growth rates lead to such high 
enzyme dilution rates that the enzymes of the "less preferred" substrates are 
diluted to near-zero levels. On the other hand, substrates that support low 
growth rates fail to efficiently dilute the enzymes of the other substrates. A 
more precise statement of this argument is given in Section I3.3[ 



Now, all the experimental results described above were obtained in batch cul- 
tures. However, since the invention of the chemostat, microbial physiologists 
have acquired extensive data on microbial growth in continuous cultures. De- 
tailed analyses of this data have revealed well-defined patterns or moti f s that 



occur in diverse microb i al species growing on various substrat es (lEglil . Il995 



Harder and Dijkhuizenl . Il982l : iKovarova-Kovar and Eglil . 119981 ). The goal of 
this work is to show that the minimal model also accounts for the patterns 
observed in continuous cultures. We begin by giving a brief preview of these 
patterns in single- and mixed-substrate cultures. 

In chemostats limited by a single substrate, steady growth can be maintained 
at all dilution rates up to the critical dilution rate at which the cells wash 
out. At subcritical dilution rates, the peripheral enzyme activities of various 
cell types a nd growth-limiting substrates invariably show one of the following 
two trends (IHarder and Dijkhuizenl . Il982l ) : 



(1) If the peripheral enzyme is fully constitutive, its activity decreases mono- 
tonically with the dilution rate (Fig.Hh). The very same trend is observed 
even if the peripheral enzyme is inducible, but the nutrient medium lacks 
substrates that can induce the synthesis of the enzyme (Fig. [lb and 
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Figure 2. The activity of inducible enzymes passe s through a maxim um at an in 



termediate dilut ion rate, (a) Activ i ties o f PtsG (ISeeto et alJ . 12004 Fig. 1) and 



/3-galactosidase ( Silver and Mateles . 1969I . Fig. 4a) during glucose- limited growth 
of E. coli K12 and lactose-limited growth of E. coli B6, respect ively, (b) Activities 



of alcohol oxidase (lEgli and Harden . Il983l . Fig. 1) and amidase ([Clarke et al.l . Il968 



Fig. 1) during methanol-limited growth of H. polymorpha and acetamide-limited 
growth of P. aeroginosa, respectively. 

Fig. Eld). 

(2) If the peripheral enzyme is inducible and the medium contains the in- 
ducing substrate, the enzyme activity passes through a maximum at an 
intermediate dilution rate (Fig. [2]). 



Clarke and coworkers were the first to observe these trends during ace tamide- 



limit ed growth of wild-type and constitutive mutants of P. aeruginosa ( jClarke et al. 



19681 ). They hypothesized that in wild-type cells, the amidase activity exhibits 



a maximum (Fig. [2Jd) due to the balance between induction and catabolite re- 
pression. Specifically, 



at low dilution rates, catabolite repression is minimal and the rate of ami- 
dase synthesis is dependent mainly on the rate at which acetamide is pre- 
sented to the bacteria. Thus, with increase in dilution rate the amidase 
specific activity under steady state conditions increases. However, above 
D = 0.30 h~^ the growth rate has increased to the point where metabolic 
intermediates are being formed at a sufficiently high rate to cause significant 
catabolite repression. At higher dilution rates catabolite repression becomes 
dominant a nd a.t D = . 6 h~^ the enzyme concentration has decreased con- 



siderably (IClarke et al.l . Il968l . p. 232). 



In constitutive mutants, the enzyme activity decreases monotonically (Fig.[Th) 
because "they lack completely the part where, in the curves for the wild- 
type strain, . . . induction is dominant." These hypotheses have subsequently 
been invoked to rationalize the similar trends found in the wild-typ e cells 
and constitu t ive rnutants of ma ny other microbial systems (reviewed in iDeanl . 
19721 : iMatinl . Il978l : iTodal . Il98lh . We shall show below that the decline of the 
enzyme activity (at all D in constitutive mutants, and at high D in wild-type 
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Figure 3. Simultaneous and preferential growth patterns in continuous cultures, 
(a) During growth of H. polymorpha on a mixture of xylose and glycerol, both sub- 
strates are consumed at all d ilution rates up to washout (calculated from Fig. 3 of 
Brinkmann and Babel . 19921 ). Closed and open symbols show the substrate concen- 
trations during single- and mixed-substrate growth, respectively, (b) During growth 
of E. coli B on a mixture of glucose and lactose, consumption of lactose ceases at 
a diluti on rate below the ( washo ut) dilution rate at which consumption of glucose 
ceases (ISilver and Matele j . Fig. 3). 
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Figure 4. Variation of the substrate concentrations and cell density dur ing mixed 



substrate growth of H. polymorpha on a mixture of glucose and methanol (lEgli et al 



19861 . Fig. 2). The total feed concentration of glucose and methanol was 5 g L in 



all the experiments. The percentages show the mass percent of methanol in the feed. 
The critical dilution rates on pure methanol (curve labeled 100%) and pure glucose 
(curve labeled 0%) are ~0.2 and ~0.5 h~^, respectively. 



cells) is entirely due to dilution, rather than catabolite repression. 



In chemostats limited by pairs of carbon sources, two types of steady state 
profiles have been observed. Importantly, these profiles correlate with, and 
in fact, can be predicted by, the substrate consumption pattern observed in 
substrate-excess batch cultures. Indeed, pairs of substrates consumed simul- 
taneously in substrate-excess batch cultures are consumed simultaneously in 
continuous cultures at all dilution rates up to washout (Fig. [3h). In contrast, 
pairs of substrates that show diauxic growth in substrate-excess batch cultures 
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are consumed simultaneously in continuous cultures only if the dilution rate 
is sufficiently small. For instance, during growth of Escherichia coli B on a 
mixture of glucose and lactose, consumption of the "less preferred" substrate, 
lactose, declines sharply at an intermediate dilution rate, beyond which only 
glucose is consumed (Fig. [3]b) . This growth pattern has been observ ed in sev- 
eral other systems (reviewed in iHarder and Dijkhuizenl . Il976l . Il982l ). but the 
most comprehensive data was obtained in studies of the methylotrophic yeasts, 
Hans enula yolymoryha and C an dida b oidinij with mixtures of methanol + glu- 
cose ( Egli et al. . 1980 . 1982bl |al. Il986l )[^ These studies, which were performed 
with various feed concentrations of glucose and methanol, revealed several 
well-defined patterns: 



(1) The sharp decline in methanol consumption at an intermediate dilution 
rate is analogous to the phenomenon of diauxic growth in batch cultures 
inasmuch as it is triggered by a precipitous drop in the activities of the 
peripheral enzymes for methanol, the "less preferred" substrate. The tran- 
sition dilution rate was empirically defined by Egli as the dilution rate at 
which the concentration of the "less preferred substrate " achieves a suffi 
ciently high value, e.g., half of the feed concentration (lEgli et al.l . Il986 
Fig. 1). 

(2) The transition dilution rate is always higher than the critical dilution 
rate on methanol (Fig. [4^) . This was referred to as the enhanced growth 
rate ejfect, since methanol was consumed at dilution rates significantly 
higher than the critical dilution rate in cultures fed with pure methanol. 

(3) The transition dilution rate varies with the feed composition — the larger 
the fraction of methanol in the feed, the smaller the transition dilution 
rate (Fig. Sh). 



These and several other general trends, discussed later in this wor k, haye been 



observed in a wide variety of bacteria and yeasts (reviewed in lEglil . Il995 
Kovarova-Kovar and Eglil . Il998l ). 



The occurrence of the very same growth patterns in such diverse organisms 
suggests that they are driven by a universal mechanism common to all species. 
Thus, we are led to consider the minimal model, which accounts for only those 
processes — induction and growth — that occur in all microbes. In earlier 
work, computational studies of a variant of the mini mal model show ed that 



1998bD . Here, 



it captures all the growth patterns discussed above (iNarang 
we perform a rigorous bifurcation analysis of the model to obtain a complete 
classification of the growth patterns at any given dilution rate and feed con- 



centrations. The analysis reveals new features, such as threshold effects, that 



^ We shall constantly appeal to this extensive data on the growth of methylotrophic 
yeasts. A detailed exposition of the m etabolism and gene regulat ion in these organ- 



isms can be found in a recent review ( Hartner and Glieder . 20061 ) 
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Figure 5. Kinetic scheme of the minimal model (lNarangl . fl998al ). 



were missed in our earlier work. It also provides simple explanations for the 
following questions: 

(1) During single-substrate growth of wild-type cells, why does the activity 
of the peripheral enzymes pass through a maximum? 

(2) During growth on mixtures of substrates, what are conditions under 
which (a) both substrates are consumed at all dilution rates up to washout 
(Fig. [3h), and (b) consumption of one of the substrates ceases at an in- 
termediate dilution rate (Fig. [3b)? 

(3) Why does the transition dilution rate exist? Why is it always larger than 
the critical dilution rate of the "less preferred" substrate, and why does 
it vary with the feed concentrations (Fig. [4^)? 

Finally, we show that under the conditions typically used in experiments, 
the physiological (intracellular) steady states are completely determined by 
only two parameters, namely, the dilution rate and the mass fraction of the 
substrates in the feed (rather than the feed concentrations per se). In fact, we 
derive explicit expressions for the steady state substrate concentrations, cell 
density, and enzyme levels, which provide physical insight into the variation 
of these quantities with the dilution rate and feed concentrations. 



2 Theory 



Fig. [5] shows the kinetic scheme of the minimal model. Here, Si denotes the i*^ 
exogenous substrate, Ei denotes the "lumped" peripheral enzymes for Si, Xi 
denotes internalized Si, and C~ denotes all intracellular components except 
and Xi (thus, it includes precursors, free amino acids, and macromolecules). 

We assume that: 
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The concentrations of the intracellular components, denoted e^, Xj, and 
c~, are based on the dry weight of the cells (g per g dry weight of cells, 
i.e., g gdw~^). The concentrations of the exogenous substrate and cells, 
denoted and c, are based on the volume of the reactor (g and 
gdw L~^, respectively). The rates of all the processes are based on the 
dry weight of the cells (g gdw~^ h^^). We shall use the term specific rate 
to emphasize this point. 

The choice of these units implies that if the concentration of any intracel- 
lular component, Z, is z g gdw~\ then the evolution of 2; in a continuous 
culture operating at dilution rate, D, is given by 




where and r~ denote the specific rates of synthesis and degradation 
of Z, and Dz is the specific rate of efflux of Z from the chemostat. It is 
shown below that the sum, Dz + {l/c){dc/dt)z, is precisely the rate of 
dilution of Z due to growth. 

The transport and peripheral catabolism of Si is catalyzed by a unique 
system of peripheral enzymes, Ei. The specific uptake rate of Si, denoted 
rg^i, follows the modified Michaelis-Menten kinetics 

"^sA — ' s,i^i 



-^s,i ~t~ Si 

where K j is a fixed constant, i.e., inducer exclusion is negligible. 

Part of the internalized substrate, denoted Xj, is converted to C~ . The 

remainder is oxidized to CO2 in order to generate energy. 

(a) The conversion of Xj to C~ and CO2 follows first-order kinetics, i.e., 

r. — • ■ 

x,i — '^X,f^l- 

(b) The fraction of Xj converted to C~ , denoted Yi, is constant. It is 
shown below that Yi is essentially identical to the yield of biomass 
on Si. This assumption is therefore tantamount to assuming that 
the yield of biomass on a substrate is unaffected by the presence of 
another substrate in the medium. 

The internalized substrate also induces the synthesis of E^, which is as- 
sumed to controlled entirely by the initiation of transcription (i.e., mech- 
anisms such as attenuation and proteolysis are neglected), 
(a) The spec ific synthesis rate of E j follows the hyperbolic Yagil &; Yagil 



kinetics (lYagil and Yagij Il97ll ) 



where K j is a fixed constant, i.e., catabolite repression is negligible. 
In the particular case of negatively controlled operons (modulated 
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by repressors), Ve^i and ai can be expressed in terms of parameters 
associated with kinetics of RNAP-promoter and repressor-operator 
binding, respectively. In other cases, ^ must be viewed as a phe- 
nomenological description of the induction kinetics. 

(b) Enzyme degradation is a first-order process, i.e., the specific rate of 
enzyme degradation is 

(c) The enzymes are synthesized at the expense of C~ , and their degra- 
dation produces C~ . 

Given these assumptions, the mass balances yield the equations 

^ = D (Sf^i - Si) - Ts^iC, (2) 

^ = r.,-r.,-(D + l|)x., (3) 
dci , ( ^ ldc\ 

^ = (nr.,, + i^r,,^) - - rl,) - (rj, - r",) - [b + i|] c", (6) 



where s/ j denotes the concentration of Si in the feed to the chemostat. 
Following Fredrickson, we observe that eqs. (l3l)-(l5|) impli citly define the spe - 



cific growth rate and the evolution of the cell density ( IFredricksonl . 119761 ) 



Indeed, since x\ + + e\ + 62 + c =1, addition of (l3])-(l5l) yields 



/ ^ 1 dc\ dc , ^. , , 



where 



r 



J2r,,-J2il-Yi)r,, (7) 



g 

i=l i=l 



is the specific growth rate. As expected, it is the net rate of substrate accu- 
mulation in the cells. It follows from ((61) that the last term in eqs. (I3l)-(l5l) 
is precisely the dilution rate of the corresponding physiological (intracellular) 
variables. 

The equations can be simplified further because k~^, the turnover time for 
Xi, is small (on the order of seconds to minutes). Thus, Xi rapidly attains 
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quasisteady state, resulting in the equations 



dsj 

dt 
dci 

Itt 
dc 

dt 



D {Sf^i - Si) - Ts^iC, (8) 

^ti - (rg + Ci, (9) 

{Tg - D) C, (10) 



Xi ^ ' ' ^ ' ^, (11) 

Tg ^ YiTs,! + F2r.,2, (12) 

where (fTT|l -(fT2l) follow from the quasisteady state relation, ~ rs,i — rx,i- 
It is worth noting that: 

(1) The parameter, 1^, which was defined as the fraction of Xj converted to 
C~, is essentially the yield of biomass on Si. Indeed, substituting the 
relation, r^^^j ~ r^^j in ((Tj) shows that after quasisteady state has been 
attained, Yi approximates the fraction of Si that accumulates in the cell. 

(2) Although we have assumed that Yi is constant, the validity of (fT2|l does 
not rest upon this assumption. It is true even if the yields are variable 
(i.e., depend on the physiological state). 

(3) Enzyme induction is subject to positive feedback. This becomes evident 
if (flB is substituted in ([Tl): The quasisteady state induction rate, r^j, 
is an increasing function of e,. The positive feedback reflects the cyclic 
structure of enzyme synthesis that is inherent in Fig. O Ei promotes the 
synthesis of Xj, which, in turn, stimulates the induction of Ei. 

We shall appeal to these facts later. 

The equations for batch growth are obtained from the above equations by let- 
ting D = 0. In what follows, we shall be particularly concerned with the initial 
evolution of the enzymes in experiments performed with sufficiently small in- 
ocula. Under this condition, the substrate consumption for the first few hours 
is so small that the substrate concentrations remain essentially constant at 
their initial levels, So,i. In the face of this quasiconstant environment, the en- 
zyme activities move toward a quasisteady state corresponding to exponential 
(balanced) growth. This motion is given by the equations 

dci 

= rt,i - {rg + ke,i) Ci, (13) 
_ _ K.iCiSci/ {Ks^i + sp^i) ^^^^ 



YiVs,iei—^ + Y2Vs,2e2-rr^ , (15) 
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obtained by letting Si ~ so,i in eqs. (fTT il - (fT2ll . The steady state(s) of these 
equations represent the quasisteady state enzyme activities attained during 
the first exponential growth phase. 

We note finally that the cells sense their environment through the ratio, 

_ Si 

(^i = > 

which characterizes the extent to which the permease is saturated with Si. 
Not surprisingly, we shall find later on that the physiological steady states 
depend on the ratio, cxj, rather than the substrate concentrations per se. It is, 
therefore, convenient to define the ratios 

For a given cell type, cij, a^^i, af^i are increasing functions of s,, Sj^o, Sf^i, re- 
spectively, and can be treated as surrogates for the corresponding substrate 
concentrations. In what follows, we shall frequently refer to cxj, a^^i, and 
as the substrate concentration, initial substrate concentration, and feed con- 
centration, respectively. 



3 Results 



3. 1 The role of regulation 



Before analyzing the foregoing equations for batch and continuous cultures, 
we pause to consider the role of regulatory mechanisms. This is important 
because the literature places great emphasis on the regulatory mechanisms, 
whereas the model neglects them completely (K j and Ve^i are assumed to be 
constants). It is therefore relevant to ask: To what extent do the regulatory 
mechanisms control gene expression? To address this question, we shall focus 
on the well-characterized lac operon in E. coli. The analysis shows that cAMP 
activation and inducer exclusion cannot explain the strong repression observed 
in experiments. 



3.1.1 cAMP activation 



Interest in cAMP as a lac regulator began when it was discovered that the lac 
induction r ate increased ~2-fo l d upon addition of cAMP t o gluc ose-limited 



cultures (jPerlman and Pastanl . Il968l : lUllmann and Monodl . Il968l ). However 



this remained a hypothesis until Epstein et al showed that the lac induction 
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Figure 6. The relative roles of dilution and cAMP activation in lac induction, 
(a) Variation of the /3-galactosidase activity with intracellular cAMP level dur- 
ing exponent ial growth of the /ac -constitutive strain E. coli X2950 on various car- 



bon sources (lEpstein et alJ.ll975l. Fig . 2). (b,c) During batch growt h of (b) E. col i 
MC4 100 ACPTinn dKuo et al.l . hwm . Fig. 2) and (c) E. coli NC3 (I Wanner et al.l . 
19781 . Fig. 1) on various carbon sources, the steady state /3-galactosidase activity 



is inversely proportional to the specific growth rate on the carbon source. More- 
over, the /3-galactosidase activity increases at most 2-fold if 5mM cAMP is added 
to the medium (■). (d) The inverse relationship is satisfied at -D > 0.4 h~^ in 
contin uous cultures of E . coli K12 W3110 grown on glucose and glucose -|- ImM 
IPTG (IKuo et al.l . l2nn,i Fig. 1). The curves in (b)-(c) and (d) show the fits to ^ 
and (fTTll . respectively. 



rate increased with the intracellular cAMP levels (Fig. [6h). They obtained 
this relationship by measuring the /3-galactosidase activities and intracellular 
cAMP levels during exponential grov^^th of the /ac-constitutive strain E. coli 
X2950 on a wide variety of carbon sources. The use of a constitutive strain 
served to eliminate inducer exclusion: In such cells, the repressor-operator 
binding is impaired so severely that a, ~ 0, and the induction rate is unaffected 
by the inducer level (rj^j ^ The diverse carbon sources provided a way of 
varying the intracellular cAMP levels: In general, the lower th e specific growth 



rate on a substrate, the higher the intracellular cAMP level ( iBuettner et al 
1973h . 



Although Fig. [6h shows a correlation between the cAMP level and the lac 
induction rate, it does not prove that there is a causal relationship between 



14 



them. In deed, since protein d egradation is negligible in substrate-excess batch 
cultures (IMandelstaml . Il958l ) , (fT3l) yields 



r+. V ■ 

Ci ^ — - ^ — -. (16) 



Thus, the steady state enzyme level depends on the ratio of the specific induc- 
tion and growth rates. Since the substrates that yield high intracellular cAMP 
levels also support low specific growth rates, it is conceivable that Fig. [6h re- 
flects the variation of the specific growth rate, r^, rather than the specific 
induction rate, 

Eq. IHM provides a simple criterion for checking if the specific induction rate 
varied in the experiments: Evidently, K.i varies if and only if Cj is not inversely 
proportional to rg. Unfortunately, we cannot check the data in Fig. [6h with 
this criterion, since the authors did not report the specific growth rates of the 
carbon sources. However , the very same exper iments have been performed by 



others (IKuo et al.l . l2003l ; I Wanner et al.l . Il978l ) , and this data shows that the 



/?-galactosidase activity is inversely proportional to Vg (Figs. [6]b,c). It follows 
that in Figs.[6)D and[6t, the /3-galactosidase activity varies almost entirely due 
to dilution. The infiuence of cAMP on lac transcription is very modest, a fact 
that is further confirmed by the very marginal (< 2-fold) improvement of the 
/3-galactosidase activity in the presence of 5mM cAMP (Figs. Eb). The same 



is lik elv to be true of the data in Fig. [6b. Thus, the data from lEpstein et al. 



19751 , which is extensively cited to support the existence of cAMP-mediated 
lac control, forces upon us just the opposite conclusion: cAMP has virtually 
no effect on the lac induction rate. 

The data from continuous cultures also implies the absence of significant 
cAMP activation. Indeed, ([9]) implies that when lac-constitutive mutants of E. 
coli are grown in continuous cultures, the steady state /3-galactosidase activity 
at sufficiently high dilution rates is given by the relation 



e. - (17) 

This relationship also holds for wild-type cells if the medium contains satu- 
rating concentrations of the gratuitous inducer, IPTG. In both cases, cAMP 
control is significant only if Cj is not inversely proportional to D. But exper- 
iments performed with /ac-constitutive mutants (Fig. and wild-type cells 
exposed to saturating IPTG levels (Fig.[6ld) show that^^-galactosidase activity 
is inversely proportional to D for sufficiently large D\2j 



^ At small the enzyme activity is significantly lower than that predicted by (fTTl) . 
Evidently, this discrepancy cannot be due to cAMP activation. It cannot be resolved 
by enzyme degradation either: Although the data can be fitted to the equation, = 
Ve^i/{D + ke,i), the values of fce,i required for these fits are physiologically implausible 
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Catabolite repression also appears to be insignificant in cell types other than 
E. coli. The amidase activity is inversely proportional to D in the constitutive 
mutant P. aeroginosa Cll (Fig. [Th). Moreover, this declining trend is not due 
to catabolite repression, as postulated by Clarke and coworkers — it reflects 
the effect of dilution. Likewise, the activity of alcohol oxidase is inversely 
proportional to D when H. polymorpha and C. boidinii are grown in glucose- 
limited chemostats. In this case, ([9]) implies that at sufficiently large dilution 
rates 

= TT77^' (IS) 

1 + ai JJ 

which is formally similar to (fTTll . 



3.1.2 Inducer exclusion 



The data shows that inducer exclusion certainly exists. When glucose or 
glucose-6-phosphate are added to a culture growing ex ponentiallv on lactose , 
the specific uptake rate of lactose declines ~2-fold (IMcGinnis and Paigenl . 



1969I . Figs. 3,5). Likewise, in the presence of glucose, the intracellular inducer 



concentration, as measu red by the accumulation of r adioactively labeled TMG, 
also decreases ~2-fold ( IWinkler and Wilson! . 119671 . Fig. 1). But this modest 
decrease cannot account for the dramatic repression of the lac operon during 
diauxic growth: The addition of glucose-6-phosphate to a culture growing on 
lactose decreases the activity of / 3-gala ctosidase ~300-fold after only 4 gener- 
ations of growth ( Hogema et al. . 19981 . Table 1). 



We conclude that cAMP activation and inducer exclusion do infiuence the lac 
induction rate, but their effect is not commensurate with the experimentally 
observed repression. We show below that the minimal model predicts complete 
transcriptional repression, despite the absence of any regulatory mechanism. 



3.2 Single-substrate growth 



In wild-type cells, the repressor-operator binding is so tight that is large 
compared to 1. This is true even in the case of glucose (aj =5-10), which is often 



treated as a substrate consumed by constitutive enzymes (iNotley-McRobb and Ferenci 



2000l . Table 3). Under these conditions, the basal induction rate is small com- 



pared to the fully induced induction rate. Hence, at all but the smallest inducer 



(~0.5 h ^). Based on the detailed study of the ptsG operon at various dilution rates, 
the discrepancy is probably due t o enhanced expre ssion of the starvation sigma 



factor, RpoS, at low dilution rates (ISeeto et al.l . 12004 Fig. 3) 
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levels, the specific induction rate can be approximated by the expression 



and the quasisteady state induction rate, obtained by substituting (iTT]) in (jlQl) . 
is 

It is useful to make this approximation because the equations become amenable 
to rigorous analysis and yield simple physical insights. In the rest of the paper, 
we shall confine our attention to these idealized perfectly inducible systems. 

We begin by considering the growth of cultures limited by a single substrate, 
say. Si. Under this condition, S2 = 62 = 0, and the steady state values of the 
remaining variable satisfy the equations 

= ^ = D(s/,i-Si)-K,ieiaiC, (20) 

= ^ = ^-^W^ {YiVs,ieiai + ei, (21) 

at Ae,i + eiCTi 

dc 

= — = (nK.ieiai - D) c, (22) 
xi ^ . (23) 

These equations have three types of steady states: ei > 0, 62 = 0, c > 0, 
denoted -Eioi; ei > 0, 62 = 0, c = 0, denoted -Eioo; and ei = 0, 62 = 0, c = 0, 
denoted -Eqoo- The first two steady states correspond t o the persistence an d 
washout steady states of the classical Monod model ([Herbert et al.l . Il956l ). 
The third steady state, -Eqoo, is a consequence of perfect inducibility — it 
exists precisely because the specific induction rate is zero in the absence of 
the inducer. 

Fig. [7] shows the bifurcation diagram for eqs. (I20ll -(l23ll. It can be inferred from 
the following facts derived in Appendix [B1 

(1) i^ooo always exists, but it is stable if and only if 



Thus, growth cannot be sustained in a chemostat if the feed concentration 
is below the threshold level, al. Under this condition 

^e,l - ^eA^ \ S ^e,l^ \ S ^6,161, 

Ke,i + eiOi Ke,i + eicr/,1 
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Figure 7. Bifurcation diagram for single-substrate growth on Si. The blue curve 
denotes the curve of the critical dilution rate, -Dc,i, and denotes the threshold 
feed concentration. The 3-D insets depict the existence and stability of the steady 
states at any given aj^i and D. Stable and unstable steady states are represented by 
full and open circles, respectively. The dashed line shows the transition dilution rate, 
Dt^i, which plays a crucial role in mixed-substrate growth (Section l3.3|) . The figure 
is not drawn to scale: Since /ce,i ~ 0.05 h~^ and Ve^i/Ke^i ~ 1 h~^, the threshold, 
al = A;e,i/(Ve,i/-^e,i), is on the order of 0.05. 

i.e., the induction rate of Ei never exceeds the degradation rate. Hence, 
as t — s> oo, ei approaches zero, and c, cxi tend toward 0,(T/ i, respectively. 
(2) i^ioo exists if and only if ct/ i > a^, in which case it is given by the 
relations, c = 0, cxi = cr/,i, and 

■ 

It is stable if and only if D exceeds the critical dilution rate 

= FiK,ieia^,i. (24) 

Evidently, Dc,i is an increasing function of cx/ 1 that is zero when i = al. 
We show below that Dc,i is the maximum specific growth rate that can 
be sustained in the chemostat (r^ < -Dc,i)- 

At first sight, this steady state poses a paradox: There are no cells in 
the chemostat, and yet, these non-existent cells have positive enzyme 
levels. The paradox disappears once it is recognized that i^ioo reflects the 
long-run behavior of a small inoculum introduced into a sterile chemostat 
operating at a/,i > al and D > -Dc,i- As t — > oo, the specific growth rate 
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of the cells approaches the maximum level, -Dc,i- However, since Vg < -Dc,i, 



dc 
dt 



{rg -D)c< - c < 0, 



so that c — > 0, 0"! — > (T/,1, and ei approaches the positive value given 
above. Thus, in contrast to -Eqoo, the cells fail to accumulate despite 
sustained induction of the enzymes because the dilution rate is higher 
than the maximum specific growth rate that can be attained. 
(3) -Eioi exists if and only if cr/,i > al and D < -Dc,i, in which case it is given 
by the relations 

= ^ (26) 



D _{D + A;e,i) [D + FiK.iiTe,!^ 

c = Yi{sf,i-si), (28) 

all of which follow from the fact that the specific substrate uptake rate 
is proportional to D, i.e., 

D , , 

= TT- (29) 

1 

Furthermore, E\q\ is stable whenever it exists. Evidently, this is the only 
steady state that allows positive cell densities to be sustained. 
The critical dilution rate, -Dc,i, is the maximum specific growth rate that 
can be maintained at steady state. To see this, observe that at -Eioi, the 
specific growth rate equals the dilution rate. Thus, (l27l) says that specific 
growth rate increases with ai, and is maximal when a\ = cr/ i. Letting 
c"i = in (EZl) and solving for D yields (l24l) . 

Fig. (I7|) shows that at any given dilution rate and feed concentration, exactly 
one of the three steady states is stable. It is therefore straightforward to deduce 
the variation of the steady states along any path in the cx/ 1, D-plane. The 
experiments are typically performed by varying either cx/^i or D, while the 
other one is held constant. We shall confine our attention to these two cases. 



3.2.1 Variation of steady states with aj^i at fixed D 

If the feed concentration is increased at a sufficiently small fixed D (horizontal 
arrow in Fig. [7]) , steady growth is feasible only if the feed concentration lies in 
the region to the right of the blue curve, where Eiqi is stable. The variation 
of the steady states in this region is given by (l25l) - (l28l) . which imply that 
as the feed concentration increases, xi, ei, ai remain constant, and c increases 
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Figure 8. Growth of E. coli ML308 in a glucose-limited chemostat. (a) When the 
feed concentration of glucose is increased at a fixed dilution rate, the residual glu- 
cose concentration (closed symbols) is constant, and the cell density (open symbols) 



increases linearly (iLendenmannl . Il994l . Appendix B). The data for D = 0.3 h~ and 
D = 0.6 h~^ is represented by the symbols. A, A and respectively, (a) When 

the dilution rate is increased at a fixed feed concentration (100 mg L~^), the residual 
substrate concentration increases monotonical ly with D, and approaches a threshold 
at small dilution rates ( Kovarova et al. . 19961 . Fig. 2). The data was obtained at 4 
temperatures. 

linearly. Thus, any increase in the feed concentration is compensated by a pro- 
portional increase in the cell density — all other variables remain unchanged. 
We shall appeal to this fact later. 

No data is available for the inducer or enzyme levels, but experiments per- 
formed with glucose-limited cultures of E. coli ML308 show that at dilution 
rates of 0.3 and 0.6 h~\ the residual substrate concentration is indeed inde- 
pendent of the feed concentration, and the cell density does increase linearly 
with the feed concentration (Fig. [8h) . No threshold was observed in these ex- 
periments because the feed concentrations were relatively high (> 1 mg L~^). 
As we show below, the threshold concentration for glucose is ~10-50 /xg L~^. 



3.2.2 Variation of steady states with D at fixed aj^i 

If the dilution rate is increased at any fixed (T/,i < al, there is no growth 
regardless of the dilution rate at which the chemostat is operated. If cx/^i > a^, 
the variation of the steady states with D is qualitatively identical to that of the 
Monod model: The persistence steady state, -Eioi, is stable for < D < -Dc,i, 
and the washout steady state, -Eioo, is stable for D > Dc,i (vertical arrow in 
Fig. [7]). We shall confine our attention to -Eioi, since Eiqq is independent of D. 



Variation of inducer and enzyme levels It follows from (l25ll - (l26l) that 
as the dilution rate increases over the range, (0,Dc,i), the inducer level in- 
creases linearly, whereas the enzyme level passes through a maximum at 
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D = ^ (YiVs^iKe^ij fce,i- The latter is consistent with the data for systems 
with inducible enzymes (Fig. [2]). 

In contrast to the hypothesis of Clarke et al., the model explains the non- 
monotonic profile of the enzyme level as the outcome of a balance between 
induction and dilution/degradation (as opposed to catabolite repression). To 
see this, observe that the steady state enzyme level is given by the ratio of the 
induction and dilution/degradation rates, i.e., 



ei 



D + keA 



Furthermore, since Xi increases with D, so does r^^, but it saturates at suffi- 
ciently large D. At low values of D, the enzyme level increases with D because 
the induction rate increases with D. At high D, the steady state enzyme level 
inversely proportional to D because the induction rate saturates. 

The validity of the above explanation rests upon the hypothesis that the induc- 
tion rate saturates at large D. This hypothesis is supported by the data for the 
lac and ptsG operons. When batch cultures of E. coli growing exponentially 
on saturating concentrations of lactose (which are physiologically identical to 
cultures growing in a lactose-limited chemostat near the critical dilution rate) 
are expo sed to 0.4 mM IPTG, there is almost no change in the /3-galactosidase 



activity (iHogema et al.l . Il998l . Table 1). Similarly, the PtsG activity of wild- 



type and ptsG-constitutive E. coli is the same in glucose-ex cess batch cultures 



and continuous cultures operating at high dilution rates (iSeeto et al.l . 12004 
Fig. 3). Induction is therefore already saturated during growth of E. coli on 
lactose or glucose, provided the dilution rate is sufficiently large. It remains to 
be seen if this hypothesis is also valid for the other systems. However, at large 
dilution rates, the enzyme activity is more or less inversely proportional to D 
for all the systems shown Fig. [21 a result expected only if the induction rate 
saturates at large D. Thus, the decline of the enzyme activity in both consti- 
tutive mutants and wild-type cells is due to dilution, rather than catabolite 
repression. 



Variation of the substrate concentration and cell density According 
to the Monod model, the residual substrate concentration is proportional to 
the dilution rate, i. e., a^ = D/r™f\ whe re r™f^ denotes the maximum specific 
growth rate on 5*1 ( IHerbert et al.l . 119561 ). In contrast, eq. f l27ll yields 



0-1 = 0-1 + 



D 
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which implies that at small dilution rates, cJi approaches the threshold level, 
cr^, and at large dilution rates, cxi is proportional to D"^. Both these properties 
follow from the relation, i = K^iCiO"! = -D/Fi, and the steady state profiles 
of the peripheral enzyme activity. Indeed, at low dilution rates, cti is constant 
because the enzyme level increases linearly. Likewise, at high dilution rates, 
(Ji ~ D"^ because ei ~ \/D. 

The experiments provide clear evidence of threshold concentrations. Kovarova 
et al. have shown that when E. coli ML308 is grown in a glucose-limited 
chemostat, the residual glucose concentration approaches threshold levels of 
20-40 /ig/L at small dilution rates (Fig. [8]b). The authors did not identify 
the specific mechanism underlying the threshold. It may exist, for example, 
simply because the maintenance requirements preclude growth. However, we 
shall show below that in some cases, the threshold exists precisely because the 
enzymes are not induced at sufficiently low values of the substrate concentra- 
tion. 



3.3 Mixed- substrate growth 

We have shown above that in chemostats limited by a single substrate, enzyme 
synthesis is abolished if the feed concentration is below the threshold level, 
and this occurs because the induction rate of the enzyme cannot match the 
degradation rate. Our next goal is to show that in mixed-substrate cultures, 
enzyme synthesis is often abolished because the induction rate of the enzyme 
cannot match its dilution rate. To this end, we shall begin with the case of 
mixed-substrate batch cultures, where this phenomenon is particularly trans- 
parent. This analysis of batch cultures will also enable us to establish the 
formal similarity between the classification of the growth patterns in batch 
and continuous cultures. 



3.3.1 Batch cultures 

3.3.1.1 Mathematical classification of the growth patterns The 

growth of mixed-substrate batch cultures is empirically classified based on 
the substrate consumption pattern (preferential or simultaneous) during the 
first exponential growth phase. The mathematical correlate of the foregoing 
empirical classification is provided by the bifurcation diagram for the equations 



~1T = ^e,!^; — 7"^^^^^ (>^iK,ieiao,i + ^2^,2620-0,2 + Ci, (30) 




Ke,l + eiCTo,! 

620-0,2 



Ke,2 + 62^0,2 



(YlV;,ieiCro,l + 1^2^,262(70,2 + K^l) 62 



(31) 
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which describe the motion of the enzymes toward the quasisteady state cor- 
responding to the first exponential growth phase. 



We gain physical insight into the foregoing equations by observing that they 
are formally similar to the Lotka-Volterra model for two competing species 

^ = aiiVi - {buNi + buN2) iVi, (32) 
= asiVs - (&2iiVi + &22iV2) N2, (33) 

where Ni is the population density of the i-th species; aiNi is the (unrestricted) 
growth rate of the i-th species in the absence of any competion; and buNf, 
bijNiNj represent the growth rate reduction due to intra- and inter-specific 
competition, respectively. Evidently, the net induction rate, K,jejCro,i/(-K^e,j + 
eiO"o,i) — fce.iCi, and the two dilution terms in (l30ll-(l3Tll are the correlates of 
the unrestricted growth rate and the two competition terms in the Lotka- 
Volterra model. Thus, the interaction between the enzymes is competitive 
(the enzymes inhibit each other via the dilution terms), despite the absence 
of any transcriptional repression (the induction rate of Ei is unaffected by 
activity of Ej,j 7^ i). 

Now, the Lotka-Volterra equations are known to yield coexistence o r extinc- 



tion steady states, depending on the values of the parameters, bij (iMurrayl . 



19891 ). The formal similarity of eqs. ([SOD-dM]) and (ESD-dMl) suggests that the 
enzymes will exhibit coexistence and extinction steady states, depending on 
the values of the physiological parameters and the initial substrate concentra- 
tions. We show below that this is indeed the case. Moreover, the coexistence 
and extinction steady states are the correlates of simultaneous and preferential 
growth, respectively. 



Eqs. (!30l)-(l3Tl) have 4 types of steady states: ei = 0, 62 = 0; ei > 0, 62 = 0; 
ei = 0, 62 > 0; and ei > 0, 62 > 0, which are denoted Eqq, Eiq, Eqi, and Eu, 
respectively. Each of these steady states has a simple biological interpretation. 
If the enzyme levels reach Eqq (resp., Eu), none (resp., both) of the two 
substrates are consumed during the first exponential growth phase. If the 
enzyme levels reach Eiq (resp., Eqi), only Si (resp., 5*2) is consumed during the 
first exponential growth phase. The bifurcation diagram shows the existence 
and stability of these steady states (and hence, the growth pattern) at any 
given (To,i and (To,2 (Fig. [9]). It can be inferred from the following facts derived 
in Appendix [Al 

(1) £^00 always exists (for all (To,i and cro,2)- It is stable if and only if 

i = — ke,i <0<^ ao,i < a^, i = l,2 
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Figure 9. Bifurcation diagram classifying the growth patterns of mixed-substrate 
batch cultures at various initial substrate concentrations, cro,i and cjo,2- The green 
and black lines represent the threshold concentrations for induction of Ei and E2 , 
respectively. The blue and red curves represent the locus of initial substrate con- 
centrations defined by the equations, r^^i = r*2 and rg^2 = respectively. The 
insets in each of the six regions show the dynamics at the corresponding initial sub- 
strate concentrations. The first inset shows the motion of the enzymes during the 
first exponential growth phase; the second one shows the resultant evolution of the 
substrate concentrations and cell density. 

i.e., the point, (o"o,i) 0-0,2), lies below the green line and to the left of the 
black line in Fig. [9l 
(2) EiQ exists if and only if ctq,! > cr^, in which case it is unique, and given 
by 



ei = , 62 = 0. 

2cro,i 

At this steady state, the cells consume only Si, and grow exponentially 
at the specific growth rate 

rnA = YiVsAeiaoA. 



It is stable if and only if 



rg,2 < rg,U 



(34) 
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i.e., (ctq,!, o"o,2) lies below the blue curve in Fig.[9l defined by the equation 



^^9,2 = rg,i. (35) 

Thus, £"10 is stable under two distinct sets of conditions. Either (T2,o is 
below the threshold level, (T2, in which case synthesis of E2 cannot be 
sustained because its induction rate cannot match its degradation rate. 
Alternatively, (To,2 is above the threshold level, but sustained synthesis 
of E2 remains infeasible because the dilution of E2 due to growth on Si 
is too large. Thus, r* 2 can be viewed as the highest specific growth rate 
tolerated by the induction machinery for E2: Synthesis of E2 is abolished 
whenever r^^i exceeds r*2- We shall refer to r*2 as the specific growth 
rate for extinction of E2. 

For any given (To,i > o"*, the exponential growth rate, r^^i, can be higher or 
lower than r* 2 (depending on the value of cro,2), but it cannot exceed the 
corresponding extinction growth rate, r*^. This refiects the biologically 
plausible fact that ei cannot be zero during exponential growth on 5*1. 
Eqi exists if and only if (Jo,2 > '^21 i'^ which case it is (uniquely) given by 
the expressions 



ei = 0, 62 = 



ke.2 \2 , 4Xe,2r 



9,2 



Y2Vs, 



2cro,2 

At this steady state, the cells consume only 6*2, and grow exponentially 
at the specific growth rate 

rg,2 = y2Vs,2e20-o,2- 

It is stable if and only if 

r*g,i < rg,2, (36) 

i.e., {(Jo,i,o'o,2) lies to the left of the red curve in Fig. [9l defined by the 
equation 

rli = rg,2- (37) 
Furthermore, 2 < ^"^,2 with equality being attained precisely when a2 = 

in which case 2 = ^"^,2 — 0- 
Exx exists and is unique if and only if i?io and i?oi are unstable, a condition 
satisfied precisely when 

^^9,2 > rg,x, r^ x > rg,2, 

i.e., (cro i,(To,2) lie between the blue and red curves in Fig. M At this 
steady state, both enzymes "coexist" because neither substrate supports 
a specific growth rate that is large enough to annihilate the enzymes for 
the other substrate. Hence, both substrates are consumed, and the cells 
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grow exponentially at the specific growth rate 



^^9,12 = ^iK.ieiO-o,! + ^2^,2620-0,2 



where 61,62 are the unique positive solutions of (l30l)- (i3T]l . 
(5) The blue curve and red curves are given by the equations 




(38) 



(39) 



respectively. Both equations define graphs of increasing functions passing 
through (cTijCTg), but the blue curve always lies below the red curve. 

Thus, the bifurcation diagram consists of 6 distinct regions such that exactly 
one of the four steady states is stable in each region. 

The bifurcation diagram implies that there are six possible growth patterns, 
depending on the initial substrate concentrations. Three of these growth pat- 
terns occur if (To,! < al or/and cro,i < Ug. If both (To,i and (To,2 are below their 
respective threshold concentrations, neither substrate is consumed. If only one 
of them is below its threshold level, say, (Jo,2 < CTg, only Si is consumed during 
the first exponential growth phase. However, growth is not diauxic because S2 
is never consumed. If (Jo,i > cr* and (To,2 > (^l, the model predicts the existence 
of three additional growth patterns. 

(1) If (To,!, o"o,2 lie in the region between the blue and green curves, 62 ap- 
proaches zero during the first exponential growth phase. Importantly, 
since cro,2 is higher than the threshold level, E2 is synthesized upon ex- 
haustion of 5*1 at the end of the first exponential growth phase. Hence, 
growth is diauxic with preferential consumption of Si. 

(2) If (To^i, (7q 2 lie in the region between the red and black curves, 61 ap- 
proaches zero during the first exponential growth phase, i.e., there is 
diauxic growth with preferential consumption of 6*2. 

(3) If cTo,i, (70,2 lie between the red and blue curves, 61 and 62 attain posi- 
tive values during the first exponential growth phase. Consequently, both 
substrates are consumed until one of them is exhausted, and there is no 
diauxic lag before the remaining substrate is consumed. 

We show below that the foregoing classification is consistent with the data. 

3.3.1.2 Growth patterns at saturating substrate concentrations 

Physiological experiments, which are generally performed with saturating sub- 
strate concentrations (o"o,i,cro,2 ~ 1), show that depending on the cell type. 
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Figure 10. Bifurcation diagrams corresponding to diauxic growth under substrate-ex- 
cess conditions (cro,i,cro,2 ~ !)• (a) Preferential consumption of Si. (b) Preferential 
consumption of 82- 



the two substrates are consumed sequentially or simultaneously. The behav- 
ior of the model is consistent with this observation. To see this, observe that 
depending on the values of the physiological parameters (which, in effect, de- 
termine the cell type), there are three possible arrangements of the red and 
blue curves, each of which yields a distinct growth pattern at saturating sub- 
strate concentrations. Indeed, 5*1 is consumed preferentially if (1, 1) lies below 
the blue curve (Fig.fTOk). This occurs precisely when the physiological param- 
eters satisfy the condition 

Si>Si+ ' fSiV (40) 

obtained from (l38l) by neglecting enzyme degradation {ke^i = 0). Likewise, S2 
is consumed preferentially if (1, 1) lies above the red curve (Fig. [TOb). i.e., the 
physiological parameters satisfy the condition 

^^''>^^ '.. feV, (41) 



Ke,2 Ke,l Y2Vs,2Ke,l2 \K,^l, 

obtained from ( l39ll upon setting ke^i = 0. Finally, the substrates are consumed 
simultaneously if (1,1) lies above the blue curve and below the red curve 
(Fig. [9]), i.e., both the foregoing conditions are violated. 



The bi ological meaning of these conditions was discussed in earlier work (iNarang and Pilyugin 



2QQ7al ). Indeed, one can check that (i40ll and f l4Ti l are equivalent to the condi- 
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tions 



a < a* = — '■ 

2 

a> a. = — y ■ 



respectively, where a = -y/l2K,2K,2/Y is a measure of the maximum 

specific growth rate on 5*2 relative to that on Si, and Ki = Ke^i/ \JVe^i/ {YiVs^i) 
is a measure of saturation constant for induction of Ei. Thus, the substrates 
are consumed preferentially whenever the maximum specific growth rate on 
one of the substrates is sufficiently large compared to that on the other sub- 
strate. Just how large this ratio must be depends on the saturation constants, 
Kj. Enzymes with small saturation constants are quasi-constitutive — their 
synthesis cannot be abolished even if the other substrate supports a large 
specific growth rate. In contrast, synthesis of enzymes with large saturation 
constants can be abolished even if the other substrate supports a comparable 
specific growth rate. 



3.3.1.3 Transitions triggered by changes in substrate concentra- 
tions Figs. [QHTOI imply that the growth pattern can be changed by altering 
the initial substrate concentrations. We show below that this conclusion is 
consistent with the data. 

As a first example, consider a cell type that prefers to consume 5*1 under 
substrate-excess conditions (cro,i,cro,2 ~ !)• Its growth pattern at any (To,i, cto,2 
is then described by Fig. [TOk . which implies that if (To,i is decreased suffi- 
ciently from 1, both substrates will be consumed simultaneously. Such behav- 
ior has been observed in experiments. If the initial concentrations of glucose 
and galactose are 20 mg/L and 5 mg/L, respectively, E. coli ML308 consumes 
glucose before galactose (Fig.fTTk). If the initial concentration of glucose is re- 
duced to 5 mg/L or lower, glucose and galactose are consumed simultaneously 
(Fig.[IIb). 

The data in Figs. [TTk.b are open to question because the initial cell densities 
are so large (~ 10^ cells L~^) that the substrates may have been exhausted 
before the enzymes reached the quasisteady state corresponding to balanced 
growth. However, similar results have also been obtained in studies performed 
with very small initial cell densities (10^-10^ cells L~^). Indeed, Schmidt and 
Alexander observed that "simultaneous use of two substrates is concentration 
dependent . . . Pseudomonas sp. ANL 50 mineralized glucose and aniline simul- 
taneously when__2resen^_aOjigL~J__b^ metabolized them diauxically at 300 



/ig L (iSchmidt and Alexanded . 119851 . Figs. 3-4). Likewise, Pseudomonas 
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Figure 11. Transitions of the growth patterns triggered by changes in the initial 
substrate concentrations. Upper panel: Batch growth of Escherichia coli ML308 
on mixtures of glucose and galactose, (a) At hi gh initial glucose co ncentrations (20 
mg L~^), glucose is consumed before galactose ( Lendenmann . 19941 ). (b) At low ini- 
tial gluc ose concentration s (5 mg L~^), glucose and galactose are consumed simulta- 
neously (lEgli et all . Il99,i Fig. 1). Lower pane l: Batch growth of Pseudomonas on 
PNP and PNP + glucose (jSchmidt et al.l . ll987l . Figs. 1,4). DPM denotes disintegra- 
tions per min of ^^C-labelled PNP and glucose, (c) No PNP is consumed if its initial 
concentration is 10 fig L~^. Significant consumption occurs at initial concentrations 
of 50 and 100 L~^. (d) PNP and glucose are consumed simultaneously if their 
initial concentrations are high (3 and 10 mg L~^, respectively). 



acidovorans consumed acetate before phenol when the concentrations of these 
substrates were >70 and 2 /ig L~^, respectively. If the initial concentration of 
acetate was reduced to 13 fig L~\ both substrates were consumed simultane- 
ously. 

In terms of the model, these transitions from sequential to simultaneous sub- 
strate consumption can be understood as follows. Sequential substrate con- 
sumption occurs at high concentrations of the preferred substrate because it 
supports such a large specific growth rate that the enzymes associated with 
the secondary substrate are diluted to extinction. Reducing the initial con- 
centration of the preferred substrate decreases its ability to support growth, 
and thus dilute the secondary substrate enzymes to extinction. Consequently, 
both substrates are consumed. 
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The growth of Pseudomonas sp. on p-nitrophenol (PNP) and glucose furnishes 
additional examples of growth pattern transitions driven by changes in the ini- 
tial substrate concentrations. Schmidt et al. found that no PNP was consumed 
if its initial concentration was 10 /xg L~^, and significant mineralization oc- 
curred at the higher initial concentrations of 50 and 100 fig (Fig. [Tib). 
It follows that the threshold concentration for PNP lies between 10 and 50 
fig L~^. Importantly, the addition of 20 mg glucose to a culture containing 
10 fig PNP failed to stimulate PNP consumption. The authors concluded 
that "the threshold was not a result of the fact that the concentration of PNP 
was too low to meet maintenance energy requirements of the organism, but 
rather supports the hypothesis that the concentration was too low to induce 
degradative enzymes." Finally, it was observed that when the initial concen- 
trations of PNP and glucose were increased to saturating levels of 3 mg 
and 10 mg respectively, both substrates were consumed simultaneously 
(Fig.fTTH). All these results are consistent with the bifurcation diagram shown 
in Fig. B 



Figs. [9]-[T0l extend the results obtained in iNarang and Pilyuginl . l2007al . There, 
we constructed bifurcation diagrams describing the variation of the growth 
patterns in response to changes in the physiological parameters at (fixed) satu- 
rating substrate concentrations. These genetic bifurcation diagrams explained 
the phenotypes of many different mutants and recombinant cells. Figs. lOHTOl 
can be viewed as epigenetic bifurcation diagrams, since they describe the vari- 
ation of the growth pattern when a given cell type, characterized by fixed 
physiological parameters, is exposed to various initial substrate concentra- 
tions. 



3.3.2 Continuous cultures 

In the experimental literature, the growth pattern of mixed-substrate contin- 
uous cultures is classified based on the manner in which the two substrates 
are consumed when the dilution rate is increased at a given pair of feed con- 
centrations. Either both substrates are consumed at all dilution rates up to 
washout, or both substrates are consumed up to an intermediate dilution rate 
beyond which only one of the substrates is consumed. We show below that the 
model predicts these growth patterns. Moreover, it can quantitatively capture 
the variation of the steady state enzyme levels, substrate concentrations, and 
cell density with respect to the dilution rate and the feed concentrations. 



3.3.2.1 Mathematical classification of the growth patterns The 

mathematical correlate of the experimentally observed growth patterns cor- 
responds to the bifurcation diagram for eqs. (i8l)-f|T2l). which describes the 
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Table 1 

Necessary and sufficient conditions for existence and stability of the steady states 
for mixed-substrate growth in continuous cultures (see Appendix C for details). 



Steady 


Defining 


Biological 


Existence 


Stability 


state 


property 


meaning 


condition(s) 


condition(s) 


-'-^OOO 


ei = = 


IMo substrate 


Alwavs 


-D/ 1 < 




c = 


consumed 
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D > Dio 




c> 
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< L> < Dc,i 
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OtiIv So 

\. y x±± y J 


9 > 
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gi = 0. P9 > 0. 


Washout 


A2 > 


A 1 < A 2, 




c = 


of £^011 




I? > A,2 


-£'111 


Pl > 0,P2 > 0, 


Both 5i, S2 


A,i,A,2 >0, 


Whenever it 




c> 


consumed 


Z)< A,i,A,2,A 


exists 


-£'110 


Pl > 0,P2 > 0, 


Washout 


A,2 > A,i, 


I?> A 




c = 


of Em 


A,i > A,2 





existence and stability of the steady states in the 3-dimensional o"/ 1, cr/ 25 -D- 
space. 

Eqs. (l8l)-(fT2l) have 7 steady states, but 5 of them are "single-substrate steady 
states" (see rows 1-5 of Tabled]). Indeed, -Eqoo? -^loi, and £"100 (resp., -Eqooi 
-Eon 5 and -Eqio) are observed in chemostats limited by 5*1 (resp., 5*2). Their 
existence in mixed-substrate growth implies that even if both substrates are 
supplied to the chemostat, there are steady states at which only one of the 
substrates is consumed. We shall see below that these steady states can become 
stable under certain conditions. Here, it suffices to observe that only two of 7 
steady states, namely, -Em (ei > 0, 62 > 0, c > 0) and i^no (ei > 0, 62 > 0, 
c > 0), are uniquely associated with mixed-substrate growth. They correspond 
to simultaneous consumption of both substrates and its washout. 

The existence and stability of the steady states are completely determined 
by five special dilution rates, namely, A,j, A,i, and A (Table 1, columns 4- 
5). These dilution rates are the analogs of the special specific growth rates 
considered in the analysis of batch cultures. Indeed, the transition dilution 
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rate for Si, 

Dt, = - K,, (42) 

is the analog of the extinction growth rate, r* -. It is the maximum dilution 
rate up to which synthesis of Ei can be sustained. The critical dilution rate, 
Z)c,i, given by the relation 



is the analog of r^^j. It is the maximum dilution rate up to which growth can 
be sustained in a chemostat supplied with only Si, and satisfies the relation, 
Dei < Dt^i with equality being obtained precisely when Ui = cr*, in which 
case, Dc^i = Dt,i = (see dashed and blue lines in Fig. [7|). Finally, the critical 
dilution rate, Dc, is defined as 

where 61,62 are the unique positive steady state solutions of (l30l)-(l3T!l with 
cro,i replaced by cr/^j. It is the analog of the function, 12. Thus, r* j, j, 12, 
and their analogs, Dt i, D^^i, are defined by formally identical expressions, 
the only difference being that (7o,i is replaced by cx/ j. 

The formal similarity of the expressions for Df^i, D^^i, and r* j, j, rg^i2 
implies that for a given set of physiological parameters, the relations, o"/ j = 
0"*, Dt^i = -Dc,2 and Dt^2 = -Dc,i, define curves on the a/^i, (j/,2-plane that 
are identical to the curves on the (To,i, cro,2-plane defined by the equations, 
<^o,i = CTj*, ^"3,1 = ^^3,2 and r*2 = r^^i. It follows that the relations, aj^i = a*, 
Dt i = -Dc,2, Dt 2 = -Dc,i, determine a partitioning of the o"/ 1, cr/2-plane into 
6 regions, which is identical to the partitioning of the bifurcation diagram for 
batch cultures, the only difference being that the coordinates are o"/i,(j/2, 
rather than o"o,i,cro,2 (Fig. [T2ll. 



Each of the six regions in Fig. [12] corresponds to unique growth pattern that 
can be inferred from the existence and stability properties summarized in 
columns 4-5 of Table 1. For instance, if (o"j i,(T/2) lies in the region between 
the green and blue curves in Fig. [12], both substrates are consumed for all 
< D < Dt2, only Si is consumed for all Dt^2 < -D < -Dc,i7 and neither 
substrate is consumed because the cells are washed out for all D > -Dc,i- To 
see this, observe that in the region of interest (between the green and blue 
curves), the special dilution rates stand in the relation 

< D,,2 < A,2 <De< < 

The existence and stability conditions listed in Table [T] therefore imply that 
as D increases, 4 of the 7 steady states play no role because they are unstable 
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Figure 12. Bifurcation diagram classifying the growth patterns of mixed-substrate 
continuous cultures at various feed concentrations, aj^i and The green and black 
lines represent the threshold concentrations for induction of Ei and E2, respectively. 
The blue and red curves represent the locus of initial substrate concentrations defined 
by the equations, L'c,! = Dt^2 and -Dc,2 = Df i, respectively. The insets in each of 
the six regions show the variation of the steady states with D at the corresponding 
feed concentrations. 

whenever they exist, or do not exist at all. Indeed: 

(1) i^ooo and i^oio exist at all D > 0, but they are always unstable because 
the stability conditions for these two steady states (Dt^i,Dt^2 > and 
Dt^i < -Dc,2, respectively) are violated in the region of interest . 

(2) i^oii exists for all < D < -Dc,2, but it is unstable whenever it exists 
because the stability condition, D > D^ i, is not satisfied in the region of 
interest. 

(3) i^iio does not exist at all since one of the existence conditions, Dt^2 > Dc^i, 
is violated in the region of interest. 



It follows that only the remaining three steady states are relevant. One can 
check by appealing to Table 1 that Em is stable for < D < 1)^ 2, -E'loi is 
stable for Dt^2 < D < -Dc,i, and £^100 is stable for D > -Dc,i- Hence, both 
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substrates are consumed for all < D < Dt^2- At D = 2, 62 becomes zero. 
Only Si is consumed for all Dt^2 < < Dc^i until c becomes zero at D = -Dc,i. 
Neither substrate is consumed for D > -Dc,i- 

Similar arguments, applied to all the regions in Fig. [121 show that there are 6 
distinct growth patterns. Three of these growth patterns occur when at least 
one of the feed concentrations is below its threshold level. 

(1) If both feed concentrations are below the threshold level, -Eqoo is always 
stable, so that neither substrate is consumed at any dilution rate. 

(2) If (7/^2 < o"2 7 but cT/^i > al, EiQi is stable (i.e., only Si is consumed) for 
all < D < -Dc,i, and £'100 is stable (i.e., the cells wash out) whenever 
D > -Dc,i- This is qualitatively similar to the behavior observed during 
single-substrate growth on Si (Fig. [7j). Thus, if the feed concentration 
of 5*2 is at sub-threshold levels, it has no effect on the behavior of the 
chemostat. 

(3) If CT/i < 0"^, and (T/^2 > (^li -^011 is stable (i.e., only S2 is consumed) 
for all < D < -Dc,2, and £'010 is stable (i.e., the cells wash out) for 
D > -Dc,2- This is qualitatively similar to the behavior of the chemostat 
during single-substrate growth on S2. 

If both feed concentrations are above their respective threshold levels, there 
are three additional growth patterns. 

(1) If cT/^i and cr/,2 lie between the red and blue curves, Em (i.e., both sub- 
strates are consumed) for < D < Dc, and Euq is stable (i.e., cells wash 
out) for D > Dc. Thus, both substrates are consumed at all dilution rates 
up to washout. 

(2) If and (T/,2 lie between the green and blue curves, Em is stable (i.e., 
both substrates are consumed) for all < -D < Dt^2, Eiqi is stable (i.e., 
only Si is consumed) for all Dt^i < D < -Dc,i, and Eiqq is stable (the cells 
wash out) for all D > -Dc,i- 

(3) If cT/^i and (T/2 lie between the black and red curves, Em is stable (i.e., 
both substrates are consumed) for all < D < Dt^i, Equ is stable (i.e., 
only S2 is consumed) for all Dt^2 < -D < -Dc,2, and Eioo is stable (the cells 
wash out) for all D > -Dc,2- 

These growth patterns are the mathematical correlates of the growth patterns 
shown in Fig. [3l 

The model predicts that the growth pattern in continuous cultures can be 
changed by altering the feed concentrations. Unfortunately, there is no experi- 
mental data at subsaturating feed concentrations (cx/^j < Ks,i) due to technical 
difficulties associated with wall growth. Henceforth, we shall confine our at- 
tention to growth at saturating feed concentrations (cr/,j ~ 1). 
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We have shown above that for a given cell type, the bifurcation diagrams for 
batch and continuous cultures are formally identical — one can be generated 
from the other by a mere relabeling of the coordinate axes. This formal iden- 
tity provides a precise mathematical explanation for the empirically observed 
correlation between the batch and continuous growth patterns of a given cell 
type. To see this, suppose the cells in question consume 5*1 preferentially in 
substrate-excess batch cultures (cro,i ~ 1). The bifurcation diagram for this 
system then has the form shown in Fig. [TOk . The continuous growth patterns 
of this system are given by the very same figure, the only difference being that 
the coordinates, ctq,!, o"o,2, are replaced by cr^ i, crj2, respectively. It follows 
that if the cells are grown in continuous cultures fed with saturating substrate 
concentrations (cx/^j ~ 1), they will consume both substrates at low dilution 
rates, and only 5*1 at all high dilution rates up to washout (Fig. [3b)- A similar 
argument shows that if the cells consume both substrates in substrate-excess 
batch culture (bifurcation diagram given by a relabeled Fig. [9]) , their growth 
in continuous cultures will be such that both substrates are consumed at all 
dilution rates up to washout (Fig. [3h). 



3.3.2.2 The model predicts the observed variation of the steady 
states with D and sj^i Although the model predicts the existence of the 
observed growth patterns, it remains to determine if the variation of the steady 
states with D and sj^i is in quantitative agreement with the data. Since the 
variation of the "single-substrate steady states" with D and Sf^i is consistent 
with the data (Section 13.21) . it suffices to focus on -Em, the only growth- 
supporting steady state uniquely associated with mixed-substrate growth. 
This steady state can be observed in systems exhibiting both the simultaneous 
and the preferential growth patterns. However, since the data for simultaneous 
systems is limited, we shall focus on preferential systems. In our discussion of 
these systems, we shall assume, furthermore, that 5*1 is the "preferred" sub- 
strate. This entails no loss of generality since the equations for Si and 5*2 (and 
the associated physiological variables) are formally identical: Interchanging 
the indices does not change the form of the equations. 

Two types of experiments can be found in the literature. 

(1) Either the dilution rate was changed at fixed feed concentrations. 

(2) Alternatively, the mass fraction of the substrates in the feed, ipi = sj^i/ (s/ 1+ 
s/,2), was changed at fixed total feed concentration, Sf^t = s/^i + s/,2, and 
(sufficiently small) dilution rate. 

The model predicts that in the first case, Em is stable at all dilution rates 
below Dt 2, the transition dilution rate for the "less preferred" substrate. In the 
second case, Em is stable at all feed fractions, except those near the extreme 
values, and 1, where the model predicts threshold effects. Such threshold ef- 
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fects have been observed when the feed contained a very small fraction of one of 
the s ubstrates (IKovar et all l2002l : iNg and Pawed . Il973l : iRudolph and Gradyl . 
20021 ). Nevertheless, we shall restrict our attention to intermediate values of 
ipi at which Em is the stable steady state. 



The steady state, Em, satisfies the equations 

= D{Sf^i - Si) - Ts^iC, 








K. 



h 



(43) 
(44) 

(45) 
(46) 



These equations do not yield an analytical solution valid for all D and 
However, we show below that they can be solved for the two limiting cases of 
low and high from which the entire solution can be pieced together. 



Before considering these limiting cases, it is useful to note that eqs. (I43l) and 
(l45l) imply that 



YiT.. 



D 



Sl) + Y2 (S/,2 - S2) 

Yi{sf^i - Si) 



Yl{Sf^l - Sl) + V2(s/,2 - S2)' 



(47) 
(48) 



The first relation says that the steady state cell density is the sum of the 
cell densities derived from the two substrates. The second relation states that 
Pi, the fraction of biomass derived from Si, equals the fraction of the total 
specific growth rate, D, supported by Si. Evidently, A < 1 with equality 
being attained only during single-substrate growth. 



Eq. (l48ll immediately yields 



Y ' 



(49) 



which explains the following well-known empirical observation: At any given 
D, the mixed-substrate specific uptake rate of Si is always lowe r than the 
single-substrate specific uptake rate of Si (reviewed in lEglil . Il995l ) . This cor- 
relation is a natural consequence of the mass balances for the substrates and 
cells, together with the constant yield property. The mixed-substrate specific 
uptake rate is smaller simply because both substrates are consumed to sup- 
port the specific growth rate, D, whereas only one substrate supports the very 
same specific growth rate during single-substrate growth. 



3.3.2.3 Low dilution rates The first limiting case corresponds to di- 
lution rates so small that both substrates are at subsaturating levels, i.e.. 
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Figure 13. Upper panel: The specific substrate uptake rate increases linearly when 
the dilution rate is increased at fixed feed concentrations. Specific uptake rates of 
(a) glucose and (b) me thanol dur i ng gr owth of H. polymorpha on various mixtures of 



glucose and methanol (lEgli et al.l . 119861 . Figs. 4-5). The percentages show the percent 
methanol in the feed. The total feed concentration of glucose and methanol was 5 
g L~^. Lower panel: If the fraction of Si in the feed is increased at a fixed dilution 
rate, the specific substrate uptake rate increases linearly with (c) Gr owth of H . 



polym orpha on a mixture of glucose and methanol at D = 0.145 h (lEgli et al 



1982bl . Fig. 5). (d) Growth of E. coli K12 on various mixtures of lactose and glucose 



at D = 0.4 h ^ (calculated from Fig. 4 of Smith and Atkinson . IQBOl ). The lines 



in the figures show the specific substrate uptake rates predicted by eq. (fST]) with 
the experimentally measured single-substrate yields, Iglu = 0.55, Imet = 0.38 in 
(a)-(c), and Yqlu = ^^lac in (d). 



Si < Ks^i <^ Sf^i. Since both substrates are almost completely consumed, (l47ll- 
(l48ll become 

c^ris/,i + y2S/,2, (50) 

Evidently, /3j is completely determined by ipi, the fraction of Si in the feed. 
Moreover, Pii^ipi) is an increasing function of ipi with /3j(0) = and A(l) = 1. 
It is identical to ipi when the yields on both substrates are the same. 



Eq. (iSTll implies that ifD is increased at fixed feed concentrations, rs,i increases 
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Figure 14. Variation of the peripheral enzyme levels: (a) Activity of /3-galactosidase 
during growth of E. coli K12 on mixtures of lac tose (ImM) plus glucose 
(2mM), glucose-6-phosphate (2mM), or glycerol (4mM) (ISmith and Atkinsonl . llQSOl . 
Fig. 3). (b) Activity of gentisate-l,2-dioxygenase and glucose-6-phosphate dehy- 
drogenase (G6P pH) during growth of P. aeruginosa on salicylate + glucose at 



D = 0.02 h-i teudolph and Gradvl . I2nn2l . Fig. 3). (c) Activities of alcohol oxi- 



dase and formaldehyde dehydrogenase during growth of H. polymorpha on various 
mixtu res of glucose and methanol at D = 0.15 h~^ and sj^j = 5 g (lEgli et al 



1982bl . Figs. 7A,C). (d) Activity of /3-galactosidase during growth of E. coli K12 on 
various mixtures of l actos e and glucose at D = 0.4 h~^ (calculated from Fig 5 of 



Smith and Atkinsonl . [l98nl ). 



linearly. Likewise, if ijji, and hence, (3i, is increased at a fixed dilution rate, 
rs,i increases linearly with j3i. Both conclusions are consistent with the data. 
Fig. Us] shows that eq. ( ISTI l provides excellent fits to the data for growth of 
H. polymorpha on glucose + methanol and E. coli on glucose + lactose. The 
validity of (ISTl) at lo w dilution rates has been demons trated for many different 
systems (reviewed in iKovarova-Kovar and Eglil . Il998l ). 



The physiological variables, Xi and e,, are completely determined by the spe- 
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cific uptake rate of Si. Indeed, (l46ll and (l44ll imply that 



(52) 



(53) 



D + Y^,,,Ke,i + (3iD 



It follows from (l53l) that: 

(1) If the dilution rate is increased at fixed feed concentrations, the steady 
state enzyme level passes through a maximum. However, at every dilution 
rate, the mixed-substrate enzyme level is always lower than the single- 
substrate enzyme level. 

(2) If the fraction of Si in the feed is increased at a fixed dilution rate, the 
activity of Ei increases with [3i. The increase is linear if D is small, and 
hyperbolic if D is large. 

Both results are consistent with the data (Fig. [T4ll. 

Based on studies with batch cultures of E. coli, it is widely believed that 
glucose and glucose-6-phosphate (G6P) are strong inhibitors of /3-galactosidase 
synthesis, whereas glycerol is a weak inhibitor. It is therefore remarkable that 
the /9-galactosidase activity at any given dilution rate is the same in three 
different continuous cultures fed with lactose + glucose, lactose + G6P, and 
lactose + glycerol (Fig. [T4k). The model provides a simple explanation for 
this data. To see this, let Si and S2 denote lactose and glucose/G6P/glycerol, 
respectively. Now, (l53ll implies that any given dilution rate, the infiuence of S2 
on the /3-galactosidase activity, ei, is completely determined by the parameter. 



which is identical to ^pl, the mass fraction of lactose in the feed, because the 
yields on lactose, glucose, G6P, and glycerol are similar. It turns out that 
the feed concentrations used in the experiments were such that ipi ~ 0.45 in 
all the experiments. Consequently, the /3-galactosidase activity at any D is 
independent of the chemical identity of 5*2. 

The steady state substrate concentrations are given by the expression 



which follows immediately from (ISTl) and (l53ll . Evidently, if the dilution rate is 
increased at fixed feed concentrations, the substrate concentrations increase. 
However, comparison with (l27ll shows that at every dilution rate, the sub- 
strate concentrations during mixed-substrate growth are lower than the cor- 



/5i 




(54) 
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Figure 15. Upper panel: Variation of the glucose and galactose concentrations with 
the dilution rate durin g single- and mixed-substrate growth of E coli ML308 on glu- 
cose or/and galactose ( Lendenmannl . 1994 : Lendenmann et al. . 1996, Appendix B). 
The residual concentrations of (a) glucose and (b) galactose during mixed-substrate 
growth (a) are lower than the corresponding levels during single-substrate growth 
(a). Lower panel: The relationship between the substrate concentration and the 
peripheral enzyme level, (a) Growth of C. boidinii on various r nixtu res of methanol 
and glucose at D = 0.14 h"^ and sj^t = 5.0 g L"^ JSgh et al.l . Il99.i Fig. 7). When 



the activity of alcohol oxidase increases (0 < Pmet ^ 0.4), the methanol concentra- 
tion is constant. When the activity of alcohol oxidase saturates (/3met ^ 0.4), the 
methanol concentration increases, (b) Residual concentrations of galactose (open 
symbols) and glucose (closed symbols) during growth on E. coli ML308 on various 
mixtures of galac tose and glucose at D = 0.3 h~^ and Sff = 1 ( A, A), 10 (0, ♦), and 
100 (□) mg ( Lendenmannl . 19941 : Lendenmann et al. . 19961 ). 



responding levels during single-substrate growth. This agrees with the data 
(Figs. fTSk.b). Eq. f l54ll also implies that if the fraction of Si in the feed, and 
hence, is increased at a fixed dilution rate, the substrate concentration 
is constant when /5j is small, and increases when /3j is large. This occurs be- 
cause at small (3i, the linear increase of j is driven entirely by the linear 
increase of ej. Since the enzyme level saturates at sufficiently large Pi, further 
improvement of r^^j is obtained by a corresponding increase of the substrate 
concentration. FigfTSb shows that this conclusion is consistent with the data. 
To be sure, there are instances in which the substrate concentrations appear 
to increase for all < /3j < 1 (Fig. fTSH). The model implies that this occurs 
because YiVs^iKe^i is so small that the peripheral enzymes are saturated at 
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relatively small values of (3iD. 



In general, one expects the mixed-substrate steady states to depend on three 
independent parameters, namely, D, a/^i, and However, the model im- 
plies that at sufficiently small dilution rates (such that Sj -C s/ j), the steady 
state specific uptake rates, and hence, the steady state enzyme levels and 
substrate concentrations, are completely determined by only two indepen- 
dent parameters, D and Experiments provide direct evidence support- 
ing this conclusion. Indeed, the data in Fig. [TSl d were obtained by perform- 
ing experiments with three different values of the total feed concentration 
{sf^t = 1, 10, 100 mg/L). It was observed that at any given fraction of galac- 
tose in the feed, the residual sugar concentrations were the same, regardless of 
the total feed concentration. Thus, the residual substrate concentrations are 
completely determined by the dilution rate and the fraction of the substrates 
in the feed (rather than the absolute values of the feed concentrations). The 
same is true of the enzyme levels (Fig.[T4b). Athough three different substrates 
(glucose, G6P, glycerol) were used in the experiments, the /3-galactosidase ac- 
tivity at any given D was the same regardless of the identity of the substrate, 
since the fraction of lactose in the feed was identical (~0.45) in all three ex- 
periments. 



3.3.2.4 High dilution rates Eq. (l54l) implies that if the dilution is in- 
creased at fixed feed concentrations, the substrate concentrations increase 
monotonically. At a sufficiently high dilution rate, S2 approaches saturating 
levels (s2 ^ -^3,2)7 while Si remains at subsaturating levels (si < Kg^i). The 
dilution rate at which S2 switches from subsaturating to saturating levels, 
denoted -D^ 2, can be roughly estimated from the relation 

^ iDs,2 + ke,2)(Ds,2P2 + Y2Vs,2Ke,2), (55) 



1^2^,2 v;,2 



obtained by letting (T2 ~ 1 in eq. (I54l) . At dilution rates exceeding -0^,2, growth 
is limited by 5*1 only. Based on our earlier analysis of single-substrate growth, 
we expect that the steady state concentrations of and the physiological 
variables are completely determined by D, and the cell density increases lin- 
early with the feed concentration of 5*1 (Fig. [8]). We show below that this is 
indeed the case. 



When S2 > Ks,2, eqs. (146]) and (|44]) yield 

rs,2 ^,262 



X2 



kx,2 kx,2 



Ve,2 -^e,2 -^e,2 , j-.-, ^r-n\ 
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(d) 



Figure 16. Near the transition dilution rate, the specific substrate uptake rates are 
independent of the feed composition. Upper panel: Model simulations. The arrows 
point in the direction of increasing ^2 in the feed. The full lines show the specific 
uptake rates predicted by numerical simulations of the model. The dashed lines 
show the specific uptake rates predicted by eqs. ([5T]l and (f57ll - ([58]l . Lower panel: 
Specific upta ke rates of gluco se and methanol during mixed-substrate growth of H. 
polymorpha (jEgli et al.l . Il986l . Figs. 4-5). 



where the second relation was obtained by appealing to the approximations, 
c"2,cr/,2 ~ 1- Evidently, 

rs,2 ~ ^,262 (57) 

is a function of D. Since YiVs^i + l2'"s,2 = D, we conclude that 



and 



Xi 



kx,i 



, ei 



Xi 



D + keA KeA + Xi 



(58) 



(59) 



are also functions of D. 



Since i and 2 are completely determined by D, the curves representing 
the specific substrate uptake rates at various feed compositions must collapse 
into a single curve. This conclusion is consistent with the data for growth of H. 
polymorpha on glucose + methanol (Figs. fTGh.d). Indeed, simulations of the 
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Table 2 

Parameter values used to simulate the growth of H. polymorpha on mixtures of 
glucose (5i) and methanol (52). The yields are based on experimental measure- 
ments (Egli et al, 1986). The remaining parameters, which have not been mea- 
sured experimentally, are based on order-of-magnitude estimates obtained from well- 
characterized systems, such as the lac operon (Shoemaker et al, 2003, Appendix A) 



i 




Vs,i (h-i) 


Ks,i (g L-i) 


K^ (h-1) 


Ve,^ (h-1) 




ke,i (h-1) 


1 


0.55 


10^ 


10-2 


10^ 


9 X lO-'' 


7 X 10"^ 


0.05 


2 


0.38 


10^ 


10-2 


10^ 


3 X lO"'' 


6 X 10-"^ 


0.05 



model, performed with the parameter values in Table 2, are in good agreement 
with the data, with the exception of two minor discrepancies (Figs.fTBk.b). The 
specific uptake rates of glucose corresponding to 77%, and 62% methanol in 
the feed are discrepant for D > 0.25 h~^ (curves labeled □ and in Fig. [T6h). 
Likewise, the single-substrate specific uptake rate of methanol deviates from 
the model prediction when D > 0.15 h~^ (curve labeled A in Fig.fTBH). Under 
these exceptional conditions, the residual methanol concentrations are so high 
that toxic effects arise, thus reducing the yields on glucose and methanol. The 
observed specific uptake rates are therefore higher than the predicted rates. 



Eqs. fl56l)-(i57[) imply that at sufficiently high dilution rates, 62 and rs^2 de- 
crease with D until they become zero at D = Dt,2 (Fig. flTb). Consequently, 
Ts^i and Ci increase with D to ensure that the total specific growth rate re- 
mains equal to D (Fig.fTTb). Both conclusions are consistent with the data for 
growth of H. polymorpha on glucose + methanol. The activities of the glucose 
enzymes, hexokinase and 6-phosphogluconate dehydrogenase, increase with D 
(Fig. flTb). while the activity of the methanol enzyme, formaldehyde dehydro- 
genase, decreases with D (Fig. [TTH). 



Since rs,\ and ei are functions of D, so is 



(60) 



In fact, the only variables that depend on the feed concentrations are the cell 
density and the concentration of Si- To see this, it suffices to observe that 
eq. (Sij) yields 



D(s/,i - si 




(61) 
(62) 



where the terms in parentheses are functions of D. As expected, c increases 
linearly with s/ i because growth is limited by S\ only. Fig. [18] shows that 
numerical simulations of the model with the parameter values in Table 2 are 
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(b) 
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Figure 17. Variation of the peripheral enzyme levels for glucose and methanol. Up- 
per panel: Model simulations. The arrows point in the direction of increasing 5*2 in 
the feed. The full lines show the specific uptake rates predicted by numerical simu- 
lations of the model. The dashed lines show the enzyme levels predicted by eqs. (f53]l 
and (f56l) . ([59]) . Lower panel: Activities of (c) the glucose enzymes, hexokinase 
and 6-phoshogluconate dehydrogenase (6PGDH), and (d) the methanol enzyme, 
formaldehyde dehydrogenase, during m ixed-substrate gro wth of H. polymorpha on 
3 g glucose and 2 g methanol ( Egli et al. . 1982al . Fig. 4). 



in good agreement with the data for growth of H. polymorpha on glucose 
methanol (Fig. SI). 



3.3.2.5 The patched limiting solutions approximate the exact so- 
lution Taken together, the two approximate solutions corresponding to low 
and high dilution rates coincide with the (exact) numerical solution at all dilu- 
tion rates, except in a small neighborhood of the dilution rate at which the two 
approximate solutions intersect (compare the dashed and full lines in Figs.fTBl- 
[TSl l. This remarkable agreement between the approximate and exact solution 
obtains because as D increases, S2 changes from subsaturating to saturating 
levels so rapidly that the change is, for all practical purposes, discontinuous. 
The switching dilution rate, Ds 2, provides a good approximation to the dilu- 
tion rate at which this near-discontinuous change occurs. The model therefore 
provides explicit formulas that approximate the specific uptake rates, sub- 
strate concentrations, cell density, and enzyme levels at all dilution rates and 
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D (h-^) 



(a) (b) 

Figure 18. The model simulations are in good agreement with the data for growth 
of H. polymorpha on glucose + methanol shown in Fig. |4l The arrows point in the 
direction of increasing ^2 in the feed. The full lines show the substrate concentrations 
and cell density predicted by numerical simulations of the model. The dashed lines 
show the substrate concentrations and cell density predicted by eqs. (f5Q]l . ([54ll and 

ii-iiD. 

saturating feed concentrations. The approximate solutions corresponding to 
low and high dilution rates are valid for D < 2 and D > Ds 2, respectively. 



3.3.2.6 Egli's transition dilution rate differs from our transition 
dilution rate There is an important difference between the transition di- 
lution rates defined by Egli and us. In our model, Dt^2 is the dilution rate 
at which 62, the activity of methanol enzymes becomes zero. This transition 
dilution rate, defined by eq. (l42l) . increases with the feed concentration of 
methanol. However, at the saturating feed concentrations used in the exper- 
iments, it is essentially independent of the feed concentration. On the other 
hand, Egli defined the transition dilution rate empirically as the dilution rate 
at which the residual methanol concentration achieved a sufficiently high level 
(e.g., half the feed concentration of methanol). The mathematical correlate of 
this empirically defined transition dilution rate is the switching dilution rate, 
2 (at which S2 switches from subsaturating to saturating levels). Consistent 
with Egli's observations, Ds,2 decreases with the fraction of methanol in the 
feed. Indeed, ( [551 ) yields 

Dt2-Ds2 

132 = Y2Vs,2Ke,2- ' 



'D,^2{Ds,2 + K,2y 



where Dt,2 is essentially constant at saturating feed concentrations. It follows 
that Ds^2 is a decreasing function of (32, and, hence, ^2- Moreover, in the 
limiting cases of feeds consisting of almost pure methanol {(32 ~ 1) and pure 
glucose {(32 ~ 0), the switching dilution rate, -0^2, tends to Dc,2 and Dt^2, 
respectively. 
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(a) (b) 

Figure 19. Yields on the substrates during growth of H. polymorpha and C. boi- 
dinii on glucose + methanol, (a) The y i eld on methanol is nearly constant at al- 



most all feed concentrations (lEgli et al 

at ~0.15 h~^. (b) In H. polymorpha, the yield on g l ucose decreases 



presence of high methanol concentrations (lEgli et 
4 Discussion 



1982bl . Fig. 4). The dilution rate is fixed 
1 glucose decreases ~20% in the 



1982al . Fig. 3). 



We have shown above that a simple model accounting for only induction and 
growth captures the experimental data under a wide variety of conditions. 
This does not prove that the model is correct unless the underlying mech- 
anism is consistent with experiments. Now, the model is based on following 
3 assumptions, (a) The yield of biomass on a substrate is constant, (b) The 
induction rate follows hyperbolic kinetics, and regulatory mechanisms have a 
negligible effect on the induction rate, (c) Each substrate is transported by a 
unique system of lumped peripheral enzymes. In what follows, we discuss the 
validity of these assumptions. 



4-1 Constant yields 



Both direct and indirect evidence show that in many cases, the yield is, for 
the most part, constant. 

Indirect evidence is obtained from experiments in which the specific growth 
and substrate uptake rates were measured during mixed-substrate growth at 
various dilution rates and feed concentrations. Given these rates, eq. (fT2ll can 
be satisfied by infinitely many Yi and Y2. It is therefore striking that in all these 
experiments, (fT2l l is satisfied by the single-substrate yie lds on the two sub 



strat e s. These include grow t h of E . coli on various sugars (ILendenmann et al 



19961 

cose 



Smith and Atkinson . 1980), organic acid s ( Narang et all . 19971 ). glu- 



3-phenylpropionic acid ( Kovarova et aL . 199?!): growth of Chelahac- 
ter heintzii on glucose + nitrilotriacetic acid ( IBally et al.l . Il994l ): and growth 
of methylotrophic yeasts on methanol + glucose (Fig. [T3b ) and methanol + 
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sorbitol (lEggeling and Sahml . 



19811). 



Direct evidence is obtained from experiments in which at least one of the 
yields was measured by using radioactively labeled substrates (Fig.dHl). These 
measurements show that the yield on a substrate remains essentially equal to 
the single-substrate yield, except when t oxic effects arise, or the fra ction of 
the substrate in the feed is below ~10% (IRudolph and Gradyl . l200ll ). 



4.2 Induction kinetics 



The induction kinetics and regulatory mechanisms can influence the behavior 
of the model. Specifically, if the induction kinetics are sigmoidal, multiple sta- 
ble steady states ari se, resulting in t he manifestation of the main tenance or 
pre-induction effect (iNarangl . Il998al : iNarang and Pilyuginl . l2007bl ). Likewise, 
introduction of regulatory effects in the model will alter its quantitative pre- 
dictions (although the analysis of the data for the lac operon shows that regu- 
latory effects play a minor role). However, the essential property of the model, 
namely, the competitive (Lotka-Volterra) interaction between the enzymes, is 
unaffected by these changes — the model still predicts the existence of ex- 
tinction and c oexistence steady states obta ined with the minimal model (see 
Section 4.1 of lNarang and Pilyuginl , l2007al ) . 



4.3 Transport kinetics 



Many substra t es are imported into the cell by two (or more) transport sys- 



tems (iFerencil . Il996l . Table 1). In general, these transport systems dominate 



substrate import at different substrate concentrations. For instance, galactose 
is transported by the low-affinity galactose permease under substrate-excess 
conditions, and the high-affinity methylga l actosi dase system under substrate- 
limiting conditions ( iWeickert and Adhyal . Il993l ). If these transport systems 
dominate in distinct regimes of the dilution rate, they are unlikely to have 
a significant effect on the model. However, if the two regimes overlap, it is 
difficult to predict the outcome without analyzing a suitably extended model. 
This will be the focus of future work. 



Lumping of peripheral enzymes has a profound effect if the inducer is syn 
thesized and consumed by coordinately expressed enzymes. Following Sav- 



ageau, we illustrate this point by appealing to the lac operon (jSavagead . 12001 



Fig. 10). It turns out that the concentration of the inducer, allolactose, is com- 
pletely determined by the concentration of intracellular lactose. Moreover, the 
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dynamics of intracellular lactose, denoted x, are given by the equation 



dx s 

— = Vien 

dt "^Ks + s 



hrj^Ci)X 



TgX, 



where s denotes extracellular lactose, denotes permease, and e?, denotes 
/?-galactosidase. Since permease and /3-galactosidase are encoded by the same 
operon, intuition suggests that they must be coordinately expressed, i.e., 
Cp = eb- It follows that the quasisteady state concentration of intracellular 
lactose, X ~ {Vs/kx)s/{Ks + s), is independent of the enzyme level. Thus, pos- 
itive feedback, a key component of the model disappears if the inducer (or its 
precursor) is synthesized by coordinately expressed enzymes. Under this con- 
dition, the positive feedback generated by one of the enzymes (e.g., permease) 
is neutralized by the other enzyme (e.g., /5-galactosidase). 

The relative rates of permease and /3-galactosidase expression are not known, 
but there is good evidence that galactosidase and transacetylase are not ex- 
pressed coordinately (reviewed in lAdhyal . l2003l ). Indeed, Ullmann and cowork- 
ers have shown that the higher the p rotein synthesis rate, the higher the rela- 



tive rate of transacetylase synthesis (jPanchin et al.l . Il98ll ). Furthermore, this 



uncoupling occurs at the level of transcription, and involves the transcription 
termination factor, Rho f lGuidi-Rontani et al.l . 11984 ). It seems likely that per- 
mease and /3-galactosidase are also uncoupled, but this remains a hypothesis 
until direct evidence is obtained. 



5 Conclusions 



We formulated a minimal model of mixed-substrate growth that accounts for 
only induction and growth (but no regulation), the two processes that occur 
in all microbes including bacteria and yeasts. 

(1) Quantitative analysis of the data for both bacteria and yeasts shows that 
the decline of the enzyme activities in constitutive and wild-type cells, 
typically attributed to specific regulatory mechanisms, is almost entirely 
due to dilution. 

(2) We analyzed the model by constructing the bifurcation diagrams for both 
batch and continuous cultures. 

(a) The bifurcation diagram for batch cultures describes the substrate 
consumption patterns at any given initial substrate concentrations. 
It shows that: (i) When the initial concentrations of one or more 
substrates are sufficiently small, threshold effects arise because the 
enzymes cannot be induced, (ii) Even if the initial concentrations of 
both substrates are at supra-threshold levels, induction of the en- 
zymes for one of the substrates may not occur because of dilution. 



48 



Thus, diauxic growth is feasible even if there are no regulatory mech- 
anisms, (c) The substrate consumption pattern can be switched from 
sequential to simultaneous (or vice versa) by altering the initial sub- 
strate concentrations. All these conclusions are consistent with the 
data in the literature, 
(b) The bifurcation diagram for continuous cultures, which describes the 
substrate consumption patterns at any given dilution rate and feed 
concentrations, is formally identical to the bifurcation diagram for 
batch cultures. This gives a precise mathematical basis for the well- 
known empirical correlation between the growth patterns in batch 
and continuous cultures. Specifically, substrates that are simultane- 
ously consumed in substrate-excess batch cultures are consumed si- 
multaneously in continuous cultures at all dilution rates up to washout. 
However, substrates that are sequentially consumed in substrate- 
excess batch cultures are consumed simultaneously in continuous cul- 
tures only if the dilution rate is sufficiently small — at large dilution 
rates, only the preferred substrate is consumed. This switch occurs 
precisely at the so-called transition dilution rate because synthesis of 
the secondary substrate enzymes is abolished by dilution. The tran- 
sition dilution rate is significantly higher the maximum growth rate 
on the secondary substrate. This "enhanced growth effect" occurs be- 
cause the enzyme levels are positive at the critical dilution rate of 
the secondary substrate — they become vanishingly small only if the 
cells are subjected to even higher dilution rates. 
(3) At the saturating feed concentrations typically used in experiments, the 
physiological steady states are completely determined by only two inde- 
pendent parameters, the dilution rate and the mass fraction of the sub- 
strate in the feed. We derived analytical expressions for the steady state 
values of the substrates, cell density, and enzymes that provide good fits 
to the extensive data for the growth of methylotrophic yeasts on glucose 
and methanol. 
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A Bifurcation diagram for batch cultures 



A.l Conditions for existence and stability of the steady states 



The steady states that can be attained during the first few hours of batch 
growth are given by the equations 

= ^ = Ve,i-^r^^ {YiVs,,e^ao,i + ^2^,262(70,2 + Ki) e,, (A.l) 

at Ae,l + eiCTo,! 

= ^ = Ve,2-^-^^ (^1^,161^0,1 + ^2^,262^0,2 + ^,2) 62- (A.2) 

at Ae^2 + 62(70,2 
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The Jacobian has the form 



^(61,62) 



Jll -1^2^,261^0,2 
-FiK,ie20-o,i J22 



(A.3) 



where 



Jll = V"e,l 



Ke.l(^0, 



J22 = V"e,2 



Ke,i + eido.i 

Ke,2<^0,2 



620-0,2 



2 - 2YiV^^,iei(To,i - ^2^,2620-0,2 - K,i, (A. 4) 
2 - ^iK,i6iO-o,i - 21^2^,2620-0,2 - K,2- (A. 5) 



In what follows, we derive the existence and stability conditions for each of 
the 4 steady states. 



(1) EoQ (ei = 0,62 = 0): It is evident from (lATTI l- llA^ll that this steady 
always exists. Since 



Ve.lO"0, 



J{Eoo) = 
Eqq is stable if and only if 



- - Ki 



V"e,2g'0,2 



-h 



e,2 



ke,i <0 ^ o-o,i < a* 



for i = 1,2. 

(2) EiQ (ci > 0, 62 = 0): This steady state exists if and only if the equation 



YiVsAeiaoA + k, 



e,l 



" Ke,l + eiO-0,1 

has a positive solution, ei > 0. Such a solution exists precisely when 

rl,i > ^ ao,i > al, 

in which case 



(A.6) 



61 



2a( 



0,1 



= FiK,i '-^ , (A.7) 
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where Vg^i = 1^1^,161(70,1 1^^^ is the exponential growth rate at Eiq. 
To determine the stability conditions, observe that 



Ju -F2V;,2eicro,2 
J22 



where 



Ju — VeA 



eio-0,1 



2FiK,ieicro,i - keA, 



J22 — V"e,2 



(^0,2 



K. 



e,2 



FiK,ieiao,i - ke,2- 



Now, Jii is always negative because (IA.6I) implies that at E 



10, 



K,i — 



KeA 



(FiK,iei(To,i + k 



e,l) , 



so that 



^11 = - ('5^iK,ieicro,i + ke,i] 



eiaoA 



KeA + ei(7o,: 



>'iK,ieioro,i < 0. 



Hence, Eiq is stable if and only if 



O"0,2 



J22 < -^^ K,2-^ < yiVs,ieiao^i\j^^^ + ke,2, 

^ r*g,2 < rg,l- 



(3) Eqi (ei = 0, 62 > 0): Arguments similar to those used above for £^10 show 
that Eqi exists if and only if 



r* > ^ ao,2 > (yl, 



in which case 



62 



- (^e,2 + 



Y2Vs, 



Ke.2 + 



Y2Vs, 



Eol 



2o-o,; 



„max _ v T/ 



YiVs, 



Y2Vs, 



' g,2 2 

where rg 2 = ^2^,262(^0,2 1^;^^^ is the exponential growth rate at -Eoi- 



(A.8) 
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(4) Ell (ei > 0,62 > 0): This steady state exists if and only if the equations 



= Ve,i - """'^ {YiVs,ieiao,i + ^2^,262^0,2 + , (A.9) 

Ke,l + CiCTo,! 

= Ve,2^^^ iYiVs,ieiao,i + ^2^,262^0,2 + ^,2) , (A.IO) 

^e,2 + e2cro,2 

have positive solutions, 61,62 > 0. The conditions for such solutions can 
be determined by eliminating one of the variables, 61, 62, from the above 
equations. To this end, observe that eqs. (IA.9l) - (|A.10p imply the relation 

Ve,l^^^ ke,l = V,,2^^^ ke,2. (A.ll) 

Ke,l + 6iCTo,i Ke,2 + 62(10,2 

Hence, we can eliminate 62 from (IA.9[ l- (IA.10p by solving f lA.llI) for 62, 
and substituting it in ( IA.9I1 to obtain the equation 

■ ke,i = yiV/,,1610-0,1 + 1214,277 ?n ; — 

Ae,i + 610-0,1 Ve^ij^^j^^ - keA + ke,2 

- Y2Vs,2Ke,2, 

which has a (unique) positive solution, 61 > 0, if and only if 

'~ ~ — kp 1 > YoVc 2~7T7 ' — YoVs, ^Kp 2- 

Ke,l ^'^ ^'^ - fce,l) + kp,2 

One can check that this is equivalent to the condition 

^1,1 = - ke,l > Y2Vs,2e2(^0,2\E,, = ^9,2, 

i.e., Eio is unstable. 

Similarly, eliminating 61 from ( lA.9ll -f lATT0l) yields 

Ve,2^ ■ kp,2 = >2K,2e2Cro,2 + 



Ke,l + 62Cro,2 Ve,2 j^^ ^ - ^6,2 + Ke,i 

-yiK,ii?e,i, 

which has a positive solution, 62 > 0, if and only if i?io is unstable. We 
conclude that En exists and is unique if and only if both Eiq and Eqi 
(exist and) are unstable. 

To determine the stability condition, observe that flA.9P implies 

^e,i— -2 = (JaVs,iei<To,i + r2K,2e20-o,2 + ke,i] 



Ke,i + eiao,i) ^^.1 + ^1^0.1 
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Hence, 



- 2YiVs,ieiao^i - 1^2^,2620-0,2 - ke,2, 
= - (YiK,ieicro,i + ^2^,2620-0.2 + ke.i) r> ^^.^"'"^ yiK,ieio-o,i, 

' Ke,l + eiO-0,1 

and similarly, 

J22 = - (ViK.ieiO-o,! + 1^2^5,2620-0,2 + 1^2^,2620-0,2- 

Ae,2 + 620-0,2 

It follows immediately that tr J < and det J > 0. Hence, En is stable 
whenever it exists. 



A. 2 Relative magnitudes of Vg^i and r 



9,1 



The exponential growth rate on Si, rg i{ao^i), is always less than r* j(cro,j), the 
specific growth rate at which Ei becomes extinct. To see this, observe that 
(IA.7l) - flA.8P can be rewritten as 



k 



which implies that 
with equality being attained when cxj = a*, in which case, Vgi = r*^ = 0. 



A. 3 Disposition of the bifurcation curves 

The lower (blue) bifurcation curve in Fig. [9] is defined by the condition 

rg,! = r*g,2- (A.13) 

Since 



the bifurcation curve defined by (IA.13I1 satisfies the equation 
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obtained from (IA.14P by replacing i with r* 3. Substituting the definition of 
r* j into the above equation yields 



Ve,l 
KeA 



10-0,1 - 0-1, 



1 + 



+ 



keA 



1 



e,2 



10-0,2 - (^2, 



K. 



e,2 



(o"o,: 



0"n 



which defines an increasing curve on the cro,i, 0-0,2-plane passing through the 
point {o'l,a2). Moreover, since the points on this curve satisfy the relation 

(0^0,1 - 0-J > (0-0,2 - O-J , 

J^eA J^e,2 



they never go above the line 



e,2 



[0-0,2 - 0-2) . 



(A.15) 



e,2 



A similar argument shows that the upper (red) bifurcation curve in Fig. [9l 
which is defined by the equation 



ra,2 



'9,1' 



(A.16) 



satisfies the relation 

Ve.2 , 



k 



K, 



e,2 



(O"0,2 



— 0-0 



e,2 



+ 



>2K,2i^e,2 
1 



e,12 



(0-0,2 - C^i 



KeA 



(0-0,1 - 0-1 



This relation also defines an increasing curve passing through {a^, o-g), but the 
curve never goes below the line defined by (IA.15I1 because 



^e,2 / *N ^ Ve,l I *N 

0-0,2 - 0-2) > (0-0,1 - 0-J . 



K. 



e,2 



i^e,l 



Hence, the curve defined by (|A.16l) always lies above the curve defined by 
dAll). 

Finally, we note that as the two bifurcation curves approach the point, (aj", 0-2), 
they are approximated by the lines 



K.i 

Ve,2 
Ke,2 



(0-0,1 - o-*) 

(0-0,2 - 0-2) 



1 + 
1 + 



ke,l \ 


"v;,2 


YiVs^Ke,,) 


Ke,2 


ke,2 \ 




Y2Vs,2Ke,2j 


Ke,l 



(0-0,2 - 0-2) 
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which almost coincide with the line defined by (1A.15I1 because fce,i -C YiVs^iKe^i. 
Thus, for typical parameter values, the curves have the quasicusp-shaped ge- 
ometry shown in Fig. [9l 



B Bifurcation diagram for continuous cultures 



The steady states of continuous cultures satisfy the equations 



dsi 

It 

ds2 

It 

= ^ 
dt 
dei 

IF 

de2 

'dt 



D - si) - r^^ic, 
D (s/,2 - si) - rs,2C, 
(r, - D) c, 

R2 = v;,2 



eio"! 



Ke,l + eiCTi 



where r„ 



'Ke,2 + e2(T2 

YiTs,! + ^2^^3,2- The Jacobian is 
-D 



{rg + fce,i) ei, 
{rg + A;e,2)e2, 



J 



-D 



ds 





dsi 

dRi 

dsi 

dR2 
dsi 



-r^,i -c 



drs 



dei 



ds2 









dr. 



s,2 



'de2i 



C 



dr„ 



ds2 

dRi 
ds2 

dR2 
ds2 



7~, ar„ ar„ 

y aei 062 

Q dRi dRi 

dei de2 

g dR2 djh 

dei de2 



(B.l) 
(B.2) 
(B.3) 
(B.4) 

(B.5) 



(B.6) 



We begin by considering the washout steady states (c = 0), since the existence 
and stability conditions for these steady states follow immediately from the 
foregoing analysis of batch cultures. 



B.l Existence and stability of washout steady states 



At a washout steady state, c = 0, Si = s/ i,S2 = s/ 2- Such a steady state 
exists if and only if there exist ei, 62 > satisfying the equations 

= = K,i r> {YiVs,ieiaf^i + F2K,2e2(X/,2 + K,i) ei, 

dt Ae,i + eicr/,1 

de2 620"; 2 
= -TT = Ve,2^ — (yiK,ieiCr/,l + Y2Vs,2e2Crf,2 + ke,2) 62, 

at Ae,2 + e2cr /^2 
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which are formally similar to eqs. (IA.l[l - (IA.2ll . the only difference being that 
CTo.i is replaced by cr/^j. 



It follows from (1B.6I1 that the Jacobian at a washout steady steady is 



J 



-D -r,,i 

-D -r,,2 

rg-D 

dRi dRi g dRi dRi 

dsi ds2 dei de2 

dR2 dR2 g djh dR^ 

dsi ds2 dei de2 



Thus, a washout steady state is stable if and only if at this steady state, 



Tg < D, and the submatrix 



R = 



dRi dRi 

dei de2 

dR2 dR2 

de\ de2 



has negative eigenvalues. This submatrix is formally similar to the matrix 
defined by eqs. (1A.3I1 - (1A.5I1 . the only difference being that cro,i is replaced by 

Given the above results, we expect the existence and stability conditions for 
the washout steady states, -Eooo, -E'loo, -E-oio, -^iio, to be formally similar to the 
existence and stability conditions for the corresponding batch culture steady 
states, namely, Eqq, Eiq, Eqi, and En. We show below that this is indeed the 
case. 

(1) Eqqq (ei = 0, 62 = 0, c = 0): It is evident that this steady state al- 
ways exists. Since rg\^^^^^ = < D, i^ooo is stable if and only if R{Eqqq) 
has negative eigenvalues. Since i?(£'ooo) is formally similar to J{Eqo), we 
conclude that ii^ooo is stable if and only if 



DtA = 



ke^i < af^i < a* 



k 



(B.7) 



for i = 1,2. 

(2) i^^ioo (ci > 0, 62 = 0, c = 0): This steady state exists if and only if 



Ke,l + 6iCr/,i 



YiVs,ieiaf^i + k 



e,l 



(B.8) 



which is formally similar to eq. (1A.6I1 . It follows that Eiqq exists if and 
only if 

A,l > ^ (7/,i > (Ti, 
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in which case 



ei 



-Bioo 



/^e,l = yiK,l ^ "^^^-'^^ ^\ ^^^-^^ (B.9) 

where Dc,i = ^sUioq ~ 1^1*^1 1 £100 ^'^^ critical dilution rate during 
single-substrate growth on Si. 

EiQQ is stable if and only ii D > Dc,i and Dt^2 < -Dc,i, where the latter 
relation ensures that the submatrix, R{Eiqo), which is formally similar to 
J(£'io), has negative eigenvalues. In the particular case of single-substrate 
growth on Si, the second stability condition is always satisfied (-Dt,2 = 
—ke,2 < < -Dc,i), SO that Eiqq is stable whenever D > -Dc,i- 
-£'010 (ci = 0, 62 > 0, c = 0): Arguments similar to those used above for 
EiQQ show that i?oio exists if and only if 

A,2 > ^ (T/,2 > (72, 

in which case 



-(^e,2 + ,|f,) + V(^e,2 + fe) + 

62 



2 , 4K^.2Dt.2 



Y2Vs,2 



lEoio 



■ 2 4K,oDt,2 



I^c,2 = Y2Vs,2— """^ ^\ (B.IO) 

where -Dc,2 = ''^glEoio ~ ^'^S''^^'^^'^\eqxq critical dilution rate during 

single-substrate growth on 5*2. It is stable if and only if D > ^£,2 and 

A,l < /^c,2. 

EiiQ (ei = 0, 62 > 0, c > 0): This steady state exists if and only if the 
equations 

= K,! ^ (VlKaeiCT/,! + ^2^,262^^,2 + K^) , (B.ll) 

Ae,i + ei(T/,i 

= Ve,2^^^ {YiVs,ieiaf,i + Y2Vs,2e2(Tf,2 + ^,2) , (B.12) 

have positive solutions. These equations are formally similar to eqs. (1A.9I1 - 
(lA.lOl) . Hence, £"110 exists and is unique if and only if 

< A,2 < D,,i, < A,i < D,^2. 

Since R^Ehq) is formally similar to J {En), the eigenvalues of R^Euq) 
are always negative. It follows that Euq is stable if and only if 

D > Dc = r-gl^^^^ = YiVs^iCiai + F2K,2e2cr2|^^^p . 
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We cannot solve explicitly for Dc, but its properties can be inferred from 
an implicit representation, which can be derived as follows. At -Eno, the 
steady state enzyme balances imply that 



■v^ ejCT/- i = — — —- — (Di i — Dc) . 



Hence, Dc satisfies the equation 

Dc + ke,l Dc + ke,2 

We shall appeal to this equation below. 

Before proceeding to the persistence steady states, we pause to consider the 
geometry and disposition of the surfaces. 



B.1.1 Geometry and disposition of the surfaces of Dt^i{af^i) and Dc,i{crf^i) : 

The geometry and disposition of the surfaces of Dc,i and Dt i follow from 
arguments similar to those used above in the analysis of batch cultures. Indeed, 
(iRQll- dBlnll imply that 

which is formally similar to f lA.12l ). It follows that Dt^i > Dc^i, and the bifurca- 
tion curves defined by the equations, Dt^2 = -Dc,i and DtA = -Dc,2, are identical 
to the bifurcation curves in Fig. [9l the only difference being that they lie on 
the (T/^i, (Tj^2-plane (as opposed to the cto,i, cro,2-plane). 



B.1.2 Geometry and disposition of the surface 0/ Dc(o"/,i, 

The surface of Dt^i intersects the surfaces of both Dc and Dc,2 along the very 
same curves. Indeed, it follows from (IB. 131) that Dt^i intersects the surfaces of 
Dc and Dt^i along the curve 

^"2^2^6 2 . 
Dt,l + fce,2 

which is identical to the equation defining the curve, D^ i = -Dc,2- A similar 
argument shows that the surface of 2 intersects the surfaces of Dc and -Dc,i 
along the same curve. 
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Differentiation of (1B.13P yields 



-L (D+fce,l)2 (D+fce,2)2 

This immediately implies that the magnitude of Z^c (relative to the magnitudes 
of Dt,i and Dc,i) is as shown in Fig. | 



B.2 Persistence steady states 

There are three persistence steady states. 

B.2.1 Eim fei > 0, 62 = 0, 0; 

This steady state exists if and only if the equations 

= D(s/,i-si)-K,ieiCTic, (B.14) 

= ^-^^ir^ (yiK,ieiai + , (B.15) 

= YiVs^idai - D. (B.16) 

have positive solutions. These equations can be explicitly solved for Si, Ci, and 
c. Indeed, f lB.161) implies that 

D 



which can be substituted in (IB.14l) - (IB.15ll to obtain 



c = Fi(s/,i-si), (B.17) 

"1 - D + fc,,iFiV;,iire,i + /^' ^ ^^^^ 

riK^ici yiVs,iv;,i ^ ^ 

Now, Ci > for all D > 0, and c > whenever si < sj i. Hence, i^ioi exists if 
and only if (IB.19I) has a solution satisfying < ai < cr/,i. It follows from the 
graph of (Ji (Fig. IB. lb ) that such a solution exists if and only if 

> al, < D < Dc^i. 

Here, Dc^i satisfies the relation 

== J' {Dc,i + K,i) (d^,i + YiVs,iKe.i 
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Figure B.l. (a) The residual substrate concentration at £'101 satisfies the relation, 
< (Ti < (T/^i, precisely when cr/^i > al and < D < Dc,i- (b) Graphical depiction 
of the conditions for existence of -Em. 

which has the positive solution given by (1B.9I 1. 

To determine the stability of -Eioi, observe that the Jacobian is similar to the 
matrix 



9s 1 

dRi 

dsi 







dei 

dRi 
dei 







Far 



Ml 

de2 

-D 








dih 

de2 



-C 



de2 



s,l 

dRi 
ds2 



-D 



which represents the Jacobian of the new system of equations obtained if we 
make the linear coordinate change 

{si, S2, ei, 62, c} — > {si, ei, Y^si + Y1S2 + c, 62, S2} . 

It follows that i^ioi is stable if and only if dR2/de2 < and the submatrix 

dsi dei 

dRi dRi 
dsi dei 

has negative eigenvalues. One can check that thi submatrix always has negative 
eigenvalues, and 



(9i?2 _ K,2'7/,2 



de2 Ke,2 
so that Eioi is stable if and only if 

D > Dt,2 



D - Ko = Dto - D 



(B.20) 
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In the particular case of single-substrate growth on 5*1, this condition is always 
satisfied, so that i^ioi is stable whenever it exists. 



B.2.2 Eon (ei = 0, e2> 0, c> 0) 

The existence and stability conditions for i^oii can be derived by methods 
analogous to those shown above for ii^ioi- 



B.2.3 Eui (ei> 0, e2> 0, c> 0) 

This steady state exists if and only if the equations 

0=D{sf,i-si)-Vs,ieiaic (B.21) 

=D{sf^2 - S2) - ^,262^20 (B.22) 

=^<^^T^ (^'iV^.ieKTi + ^2^,262(72 + ke,l) (B.23) 

=Ve.2^^ (^K.ieiai + 1^2^,262^2 + A;e,2) (B.24) 

=YiVs,ieiai + Y2Vs,2e2a2 - D (B.25) 



have positive solutions, Sj,ej,c > 0. It follows from (IB.21l) - (IB.22l) that c > 
if and only if Si < Sf^i. 

To prove the existence conditions, it is useful to simplify the above equations. 
To this end, observe that c can be eliminated from eqs. (IB.21I )- (1B.22[ ) to obtain 
the relation 

^,2620-2 - si) = K,ieiai (s/,2 - S2) , (B.26) 
and eqs. (1B.23[ 1- (1B.25[ 1 immediately yield 

V ■ K 

' = 1,2, (B.27) 



D + ke,i (Ti 

We conclude that Em exists if and only there exist < Sj < Cj > 
satisfying (1B.25[) - (IB.27I) . We show below that such Si and Cj exist if and only 



ifO<D< A,i,A,2,/^c. 

To prove the sufficiency of the condition, < D < Dt i,Dt 2,Dc, substitute 
( IB. 271 ) in ( IB. 2511 and ( IB. 2611 to obtain the equations 

= /i(ai, (72)^^ + ^ = 1, (B.28) 

(Tl (J2 

= /2(f^i, 02) = ^ [CTI - Oil) (s/,2 - S2) - ^ {cr2 - (L2) - Si) , 

(B.29) 
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where 



1 1' Sji-^^e,! 



Now, ( 1B.28P defines a line on the di, cr2-plane with the intercepts, ai and 0^2 
(full line in Fig. IB.lb ). As D increases, it moves further from the origin, and 
passes through the point, (cx/^i, cT/,2), precisely when D = Dc. On the other 
hand, (1B.29I1 defines a curve which has the geometry of the dashed curve in 
Fig. IB.lb . To see this, observe that (011,0:2) and (o'/^i,o"/2) lie on the curve. 
Moreover, (0:1,0:2) lies to the left and below (o"/,i, o"/,2) because 

Ve i'^f i 

D < Dt,i = - ^' - ke,i oli < o-/,i, 
and the slope of the curve is positive since 

dh/da, feSf ft - + Kir - -2) 



> 



for all Qii < ai < aj^i. Finally, (0:^,0:2) lies below the growth isocline because 

£1 £2 _ YlVs,lKe,l + Y2Vs,2Ke,2 

ai W2 D + YiVs,iK,,i + Y2Vs,2Ke,2 

for all > 0. It follows that the curves defined by (IB. 281 ) and (IB. 291 ) intersect 
at a unique point, say, (cri,a2), such that o_i < ai < Of^i. Since 

V ■ K ■ 

ai<a,^ei = — — — > 

D + fee i (Ti 



Em exists and is unique. 

The necessity of the existence conditions follows from the fact that for all 

< ai < aj^i, 



V ■ K K 



D + ke,i af^i ai {D + ke 



(A, - D) . 



Hence, if Dt^i < or D > Dt^i, Em does not exist because Cj < 0. Likewise, 
if D > Dc, Em does not exist because there are no < Sj < sj^i satisfying 
(E28D. 

We conclude that Em exists if and only if 

0<D<Dt,i,Dt,2,Dc. 
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We were unable to prove the stability condition for Em. However, extensive 
numerical simulations suggest that it is stable whenever it exists. 
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Abstract 

During batch growth on mixtures of two growth-hmiting substrates, microbes con- 
sume the substrates either sequentially or simultaneously. These growth patterns 
are manifested in all types of bacteria and yeasts. The ubiquity of these growth 
patterns suggests that they are driven by a universal mechanism common to all 
microbial species. In previous work, we showed that a minimal model accounting 
only for enzyme induction and dilution explains the phenotypes observed in batch 
cultures of various wild- type and mutant /recombinant cells. Here, we examine the 
extension of the minimal model to continuous cultures. We show that: (1) Sev- 
eral enzymatic trends, usually attributed to specific regulatory mechanisms such as 
catabolite repression, are completely accounted for by dilution. (2) The bifurcation 
diagram of the minimal model for continuous cultures, which classifies the substrate 
consumption pattern at any given dilution rate and feed concentrations, provides a 
a precise explanation for the empirically observed correlation between the growth 
patterns in batch and continuous cultures. (3) Numerical simulations of the model 
are in excellent agreement with the data. The model captures the variation of the 
steady state substrate concentrations, cell densities, and enzyme levels during the 
single- and mixed-substrate growth of bacteria and yeasts at various dilution rates 
and feed concentrations. (4) This variation is well approximated by simple analyt- 
ical expressions that furnish physical insights into the steady states of continuous 
cultures. Since the minimal model describes the behavior of the cells in the absence 
of any regulatory mechanisms, it provides a framework for rigorously quantitating 
the effect of these mechanisms. We illustrate this by analyzing several data sets from 
the literature. 



Key words: Mathematical model, mixed substrate growth, substitutable 
substrates, chemostat, gene expression. 
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1 Introduction 



The mechanisms of gene regulation are of fundamental importance in biology. 
They play a crucial role in development by determining the fate of isogenic 
embroyonic cells, and there is growing belief that the diversity of biological 
organisms reflects the variation of their regulatory mechanisms (??). 

Many of the key principles of gene regulation, such as positive, negative, and 
allosteric control, were discovered by studying the growth of bacteria and 
yeasts in batch cultures containing one or more growth-limiting substrates. 
A model system that played a particularly important role is the growth of 
Escherichia coli on lactose and a mixture of lactose and glucose (?). 

It turns out that the enzymes catalyzing the transport and peripheral catabolism 
of lactose, such as lactose permease and /9-galactosidase, are synthesized or in- 
duced only if lactose is present in the environment. The molecular mechanism 
of induction was discovered by Monod and coworkers (?). They showed that 
the genes encoding the peripheral enzymes for lactose are contiguous and 
transcribed sequentially, an arrangement referred to as the lac operon. In the 
absence of lactose, the lac operon is not transcribed because the lac repressor 
is bound to a specific site on the lac operon called the operator. This prevents 
RNA polymerase from attaching to the operon and initiating transcription. In 
the presence of lactose, transcription of lac is triggered because allolactose, a 
product of /?-galactosidase, binds to the repressor, and renders it incapable of 
binding to the operator. 

When Escherichia coli is grown on a mixture of glucose and lactose, lac tran- 
scription, and hence, the consumption of lactose, is somehow suppressed until 
glucose is exhausted. During this period of preferential growth on glucose, the 
peripheral enzymes for lactose are diluted to very small levels. The consump- 
tion of lactose begins upon exhaustion of glucose, but only after a certain lag 
during which the lactose enzymes are built up to sufficiently high levels. Thus, 
the cells exhibit two exponential growth phases separated by an intermediate 
lag. Monod referred to this phenomenon as diauxic or "double" growth (?). 

Two major molecular mechanisms have been proposed to explain the repres- 
sion of lac transcription in the presence of glucose: 

(1) Inducer exclusion (?): In the presence of glucose, enzyme IIA^'^, a periph- 
eral enzyme for glucose, is dephosphorylated. The dephosphorylated IIA^'^ 
inhibits lactose uptake by binding to lactose permease. This reduces the 
intracellular concentration of allolactose, and hence, the transcription rate 
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of the lac operon. 

Genetic evidence suggests that phosphorylated IIA^^'^ activates adenylate 
cyclase, which catalyzes the synthesis of cyclic AMP (cAMP). Since glu- 
cose uptake results in dephosphorylation of IIA^^'^, one expects the cAMP 
level to decrease in the presence of glucose. This forms the basis of the 
second mechanism of lac repression. 
(2) cAMP activation (?): It has been observed that the recruitment of RNA 
polymerase to promoter is inefficient unless a protein called catabolite ac- 
tivator protein (CAP) is bound to a specific CAP site on the promoter. 
Furthermore, CAP has a low affinity for the CAP site, but when bound 
to cAMP, its affinity increases dramatically. The inhibition of lac tran- 
scription by glucose is then explained as follows. 

In the presence of lactose alone (i.e., no glucose), the cAMP level is 
high. Hence, CAP becomes cAMP-bound, attaches to the CAP site, and 
promotes transcription by recruiting RNA polymerase. When glucose is 
added to the culture, the cAMP level decreases by the mechanism de- 
scribed above. Consequently, CAP, being cAMP-free, fails to bind to the 
CAP site, and lac transcription is abolished. 

Mechanistic mathematical models of the glucose-lactose diauxie appeal to both 
these mechanisms (????) 

Over the last few decades, microbial physiologists have accumulated a vast 
amount of data showing that diauxie growth is ubiquitous (reviewed in ???). 
It occurs in diverse microbial species and numerous pairs of substitutable sub- 
strates (satisfying identical nutrient requirements). In many of these systems, 
the above mechanisms play no role. For instance, cAMP is not detectable in 
gram-positive bacteria (?), and has little effect on the transcription of pe- 
ripheral enzymes in pseudomonads (?) and yeasts (?). It has been proposed 
that in these cases, the mechanisms of repression are analogous to those of 
the lac operon, but involve different transcription factors, such as CcpA in 
gram-positive bacteria (?), Crc in pseudomonads (?), and MigA in yeasts (?). 

Paradoxically, there is growing evidence that cAMP activation and inducer 
exclusion are not sufficient for explaining lac repression. Not long after the 
discovery of cAMP activation (??), Ullmann presented several lines of evi- 
dence showing that this mechanism played a relatively modest role in diauxie 
growth (?). But the most compelling evidence was obtained by Aiba and 
coworkers, who showed that the cAMP levels are essentially the same during 
the first and second growth phases; furthermore, glucose-mediated lac repres- 
sion persists in the presence of high (5mM) exogenous cAMP levels, and in 
mutants of E. coli in which the ability of cAMP to influence lac transcription 
is completely abolished (??). In Section ??, we shall show that the positive 
correlation between lac expression and intracellular cAMP levels (?), which 
forms the foundation of the cAMP activation mechanism, does not reflect a 
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causal relationship: The observed variation of lac expression is almost entirely 
due to dilution (rather than cAMP activation). 

The persistence of lac repression in cAMP-independent cells has led to the 
hypothesis that inducer exclusion is the sole cause of repression (?). How- 
ever, this mechanism exerts a relatively mild effect. The lactose uptake rate is 
inhibited only ~50% in the presence of glucose (?). 

Thus, transcriptional repression, by itself, cannot explain the several hundred- 
fold repression of lac during the first exponential growth phase of diauxic 
growth. In earlier work, we have shown that complete repression is predicted 
by a minimal model accounting for only induction and growth (??). 

Yet another phenomenon that warrants an explanation is the empirically ob- 
served correlation between the substrate consumption pattern and the spe- 
cific growth on the individual substrates. For instance, in the case of diauxic 
growth, it has been observed that: 

In most cases, although not invariably, the presence of a substrate per- 
mitting a higher growth rate prevents the utilization of a second, 'poorer' 
substrate in batch culture (?). 

It turns out that sequential consumption of substrates is not the sole growth 
pattern in batch cultures. Monod observed that in several cases of mixed- 
substrate growth, there was no diauxic lag (??), and subsequent studies have 
shown that both substrates are often consumed simultaneously (?). The oc- 
currence of simultaneous substrate consumption also appears to correlate with 
the specific growth rates on the individual substrates. Based on a comprehen- 
sive review of the literature, Egli notes that: 

Especially combinations of substrates that support medium or low maxi- 
mum specific growth rates are utilized simultaneously (?). 

Recently, we have shown that the minimal model provides a natural explana- 
tion for the foregoing correlations, which hinges upon the fact that the enzyme 
dilution rate is proportional to the specific growth rate (?). Roughly, substrates 
that support high growth rates lead to such high enzyme dilution rates that 
the enzymes of the "less preferred" substrates are diluted to near-zero levels. 
On the other hand, substrates that support low growth rates fail to efficiently 
dilute the enzymes of the other substrates. A more precise statement of this 
argument is given in Section ??. 

Now, all the experimental results described above were obtained in batch cul- 
tures. However, since the invention of the chemostat, microbial physiologists 
have acquired extensive data on microbial growth in continuous cultures. De- 
tailed analyses of this data have revealed well-defined patterns or motifs that 
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Figure 1. If the induction rate of the peripheral enzyme is independent of the 
inducer level, its activity decreases with the dilution rate, (a) The activities of 
/3-galactosidase (?, Fig. 1) and amidase (?, Fig. 2) during growth of the consti- 
tutive mutants, E. coli B6b2 and P. aeruginosa Cll, on lactose and acetamide, 
respectively, (b) The activity of formaldehyde dehydrogenase (FDH) during glu- 
cose-limited growth of H. polymorpha and C. boidinii (?, Fig. lc,d). The curves in 
(a) and (b) shows the fits to eqs. (??) and (??), respectively. 




FigTire 2. The activity of inducible enzymes passes through a maximum at an in- 
termediate dilution rate, (a) Activities of PtsG (?, Fig. 1) and /3-galactosidase (?, 
Fig. 4a) during glucose-limited growth of E. coli K12 and lactose-limited growth of 
E. coli B6, respectively, (b) Activities of alcohol oxidase (?, Fig. 1) and amidase (?, 
Fig. 1) during methanol-limited growth of H. polymorpha and acetamide-limited 
growth of P. aeroginosa, respectively. 

occur in diverse microbial species growing on various substrates (???). The 
goal of this work is to show that the minimal model also accounts for the 
patterns observed in continuous cultures. We begin by giving a brief preview 
of these patterns in single- and mixed-substrate cultures. 

In chemostats limited by a single substrate, steady growth can be maintained 
at all dilution rates up to the critical dilution rate at which the cells wash 
out. At subcritical dilution rates, the peripheral enzyme activities of various 
cell types and growth-limiting substrates invariably show one of the following 
two trends (?): 
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(1) If the peripheral enzyme is fully constitutive, its activity decreases mono- 
tonically with the dilution rate (Fig. ??a). The very same trend is ob- 
served even if the peripheral enzyme is inducible, but the nutrient medium 
lacks substrates that can induce the synthesis of the enzyme (Fig. ??b 
and Fig. ??d). 

(2) If the peripheral enzyme is inducible and the medium contains the in- 
ducing substrate, the enzyme activity passes through a maximum at an 
intermediate dilution rate (Fig. ??). 

Clarke and coworkers were the first to observe these trends during acetamide- 
limited growth of wild-type and constitutive mutants of P. aeruginosa (?). 
They hypothesized that in wild-type cells, the amidase activity exhibits a 
maximum (Fig. ??b) due to the balance between induction and catabolite 
repression. Specifically, 

at low dilution rates, catabolite repression is minimal and the rate of ami- 
dase synthesis is dependent mainly on the rate at which acetamide is pre- 
sented to the bacteria. Thus, with increase in dilution rate the amidase 
specific activity under steady state conditions increases. However, above 
D — 0.30 h~^ the growth rate has increased to the point where metabolic 
intermediates are being formed at a sufficiently high rate to cause significant 
catabolite repression. At higher dilution rates catabolite repression becomes 
dominant and at D = 0.6 h~^ the enzyme concentration has decreased con- 
siderably (?, p. 232). 

In constitutive mutants, the enzyme activity decreases monotonically (Fig. ??a) 
because "they lack completely the part where, in the curves for the wild-type 
strain, . . . induction is dominant." These hypotheses have subsequently been 
invoked to rationalize the similar trends found in the wild-type cells and con- 
stitutive mutants of many other microbial systems (reviewed in ???). We shall 
show below that the decline of the enzyme activity (at all D in constitutive 
mutants, and at high D in wild-type cells) is entirely due to dilution, rather 
than catabolite repression. 

In chemostats limited by pairs of carbon sources, two types of steady state 
profiles have been observed. Importantly, these profiles correlate with, and 
in fact, can be predicted by, the substrate consumption pattern observed in 
substrate-excess batch cultures. Indeed, pairs of substrates consumed simul- 
taneously in substrate-excess batch cultures are consumed simultaneously in 
continuous cultures at all dilution rates up to washout (Fig. ??a). In contrast, 
pairs of substrates that show diauxic growth in substrate-excess batch cul- 
tures are consumed simultaneously in continuous cultures only if the dilution 
rate is sufficiently small. For instance, during growth of Escherichia coli B 
on a mixture of glucose and lactose, consumption of the "less preferred" sub- 
strate, lactose, declines sharply at an intermediate dilution rate, beyond which 
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Figure 3. Simultaneous and preferential growth patterns in continuous cultures, 
(a) During growth of H. polymorpha on a mixture of xylose and glycerol, both 
substrates are consumed at all dilution rates up to washout (calculated from Fig. 3 
of ?). Closed and open symbols show the substrate concentrations during single- and 
mixed-substrate growth, respectively, (b) During growth of E. coli B on a mixture 
of glucose and lactose, consumption of lactose ceases at a dilution rate below the 
(washout) dilution rate at which consumption of glucose ceases (?, Fig. 3). 
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Figure 4. Variation of the substrate concentrations and cell density during mixed- 
substrate growth of H. polymorpha on a mixture of glucose and methanol (?, Fig. 2). 
The total feed concentration of glucose and methanol was 5 g L"-*^ in all the experi- 
ments. The percentages show the mass percent of methanol in the feed. The critical 
dilution rates on pure methanol (curve labeled 100%) and pure glucose (curve labeled 
0%) are ~0.2 and ~0.5 h~^, respectively. 

only glucose is consumed (Fig. ??b). This growth pattern has been observed 
in several other systems (reviewed in ??), but the most comprehensive data 
was obtained in studies of the methylotrophic yeasts, Hansenula polymorpha 
and Candida boidinii, with mixtures of methanol + glucose (????). ^ These 
studies, which were performed with various feed concentrations of glucose and 
methanol, revealed several well-defined patterns: 



^ We shall constantly appeal to this extensive data on the growth of methylotrophic 
yeasts. A detailed exposition of the metabolism and gene regulation in these organ- 
isms can be found in a recent review (?). 
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(1) The sharp decline in methanol consumption at an intermediate dilution 
rate is analogous to the phenomenon of diauxic growth in batch cultures 
inasmuch as it is triggered by a precipitous drop in the activities of the 
peripheral enzymes for methanol, the "less preferred" substrate. The tran- 
sition dilution rate was empirically defined by Egli as the dilution rate 
at which the concentration of the "less preferred substrate" achieves a 
sufficiently high value, e.g., half of the feed concentration (?, Fig. 1). 

(2) The transition dilution rate is always higher than the critical dilution 
rate on methanol (Fig. ??a). This was referred to as the enhanced growth 
rate effect, since methanol was consumed at dilution rates significantly 
higher than the critical dilution rate in cultures fed with pure methanol. 

(3) The transition dilution rate varies with the feed composition — the larger 
the fraction of methanol in the feed, the smaller the transition dilution 
rate (Fig. ??a). 

These and several other general trends, discussed later in this work, have been 
observed in a wide variety of bacteria and yeasts (reviewed in ??). 

The occurrence of the very same growth patterns in such diverse organisms 
suggests that they are driven by a universal mechanism common to all species. 
Thus, we are led to consider the minimal model, which accounts for only those 
processes — induction and growth — that occur in all microbes. In earlier 
work, computational studies of a variant of the minimal model showed that 
it captures all the growth patterns discussed above (?). Here, we perform a 
rigorous bifurcation analysis of the model to obtain a complete classification 
of the growth patterns at any given dilution rate and feed concentrations. The 
analysis reveals new features, such as threshold effects, that were missed in our 
earlier work. It also provides simple explanations for the following questions: 

(1) During single-substrate growth of wild-type cells, why does the activity 
of the peripheral enzymes pass through a maximum? 

(2) During growth on mixtures of substrates, what are conditions under 
which (a) both substrates are consumed at all dilution rates up to washout 
(Fig. ??a), and (b) consumption of one of the substrates ceases at an in- 
termediate dilution rate (Fig. ??b)? 

(3) Why does the transition dilution rate exist? Why is it always larger than 
the critical dilution rate of the "less preferred" substrate, and why does 
it vary with the feed concentrations (Fig. ??a)? 

Finally, we show that under the conditions typically used in experiments, 
the physiological (intracellular) steady states are completely determined by 
only two parameters, namely, the dilution rate and the mass fraction of the 
substrates in the feed (rather than the feed concentrations per se). In fact, we 
derive explicit expressions for the steady state substrate concentrations, cell 
density, and enzyme levels, which provide physical insight into the variation 
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Figure 5. Kinetic scheme of the minimal model (?). 



of these quantities with the dilution rate and feed concentrations. 



2 Theory 

Fig. ?? shows the kinetic scheme of the minimal model. Here, Si denotes the 
i^^ exogenous substrate, Ei denotes the "lumped" peripheral enzymes for Si, Xi 
denotes internalized Si, and C" denotes all intracellular components except Ei 
and Xi (thus, it includes precursors, free amino acids, and macromolecules). 

We assume that: 

(1) The concentrations of the intracellular components, denoted Cj, Xi, and 
c~, are based on the dry weight of the cells (g per g dry weight of cells, 
i.e., g gdw^^). The concentrations of the exogenous substrate and cells, 
denoted and c, are based on the volume of the reactor (g and 
gdw L^^ respectively). The rates of all the processes are based on the 
dry weight of the cells (g gdw~^ h~^). We shall use the term specific rate 
to emphasize this point. 

The choice of these units implies that if the concentration of any intracel- 
lular component, Z, is z g gdw~^, then the evolution of z in a continuous 
culture operating at dilution rate, D, is given by 



where r+ and r~ denote the specific rates of synthesis and degradation 
of Z, and Dz is the specific rate of efflux of Z from the chemostat. It is 
shown below that the sum, Dz + {1/ c){dc/ dt)z, is precisely the rate of 
dilution of Z due to growth. 
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(2) The transport and peripheral catabolism of Si is catalyzed by a unique 
system of peripheral enzymes, Ei. The specific uptake rate of S'^, denoted 
rs,i, follows the modified Michaelis-Menten kinetics 



rs,i — Vs,iei , 

where Vs,i is a fixed constant, i.e., inducer exclusion is negligible. 

(3) Part of the internalized substrate, denoted Xj, is converted to C~ . The 
remainder is oxidized to CO2 in order to generate energy. 

(a) The conversion of Xi to C" and CO2 follows first-order kinetics, i.e., 

r- — Ic ■ ■ 
X,l — '^X,l-'^l- 

(b) The fraction of converted to C^, denoted Kj, is constant. It is 
shown below that Yi is essentially identical to the yield of biomass 
on Si. This assumption is therefore tantamount to assuming that 
the yield of biomass on a substrate is unaffected by the presence of 
another substrate in the medium. 

(4) The internalized substrate also induces the synthesis of Ei, which is as- 
sumed to controlled entirely by the initiation of transcription (i.e., mech- 
anisms such as attenuation and proteolysis are neglected). 

(a) The specific synthesis rate of Ei follows the hyperbolic Yagil &: Yagil 
kinetics (?), 

where Ve^i is a fixed constant, i.e., catabolite repression is negligible. 
In the particular case of negatively controlled operons (modulated 
by repressors), and ttj can be expressed in terms of parameters 
associated with kinetics of RNAP-promoter and repressor-operator 
binding, respectively. In other cases, (??) must be viewed as a phe- 
nomenological description of the induction kinetics. 

(b) Enzyme degradation is a first-order process, i.e., the specific rate of 
enzyme degradation is 

(c) The enzymes are synthesized at the expense of C~, and their degra- 
dation produces C~ . 
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Given these assumptions, the mass balances yield the equations 
dsi 

= ^ isf,i - Si) - rs,iC, (2) 
^^'--'--('' + ^^)^- 

dci , ( ^ ldc\ , , 

^ = <.-r--(c + --je., (4) 

^ = {Y,r,, + Y,r.,^ - (r+ - r;,) - - r,:,) ~ [d + i|) c", (6) 

where s/^j denotes the concentration of in the feed to the chemostat. 

Following Fredrickson, we observe that eqs. (??)-(??) implicitly define the 
specific growth rate and the evolution of the cell density (?). Indeed, since 
Xi + ^2 + ei + 62 + c~ = 1, addition of (??)-(??) yields 

1 dc\ dc , 

where 

2 2 

rg = J2 - - (7) 

i=l i=l 

is the specific growth rate. As expected, it is the net rate of substrate accu- 
mulation in the cells. It follows from (??) that the last term in eqs. (??)-(??) 
is precisely the dilution rate of the corresponding physiological (intracellular) 
variables. 

The equations can be simplified further because k^}, the turnover time for 
Xi, is small (on the order of seconds to minutes). Thus, Xi rapidly attains 
quasisteady state, resulting in the equations 



dsj 

dt 
dcj 

dt 

dc 

di 



D{sf^i- Si) -rs,iC, (8) 
<i - (rg + ke,i) ei, (9) 
{rg-D)c, (10) 



^ _ Vs,ieiSi/ {Ks,i + Si) 

Xi , [LL) 

Tg ^ Yirs,i + Y2rs,2: (12) 
where (??)-(??) follow from the quasisteady state relation, ~ rs,i — rx,i- 
It is worth noting that: 
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(1) The parameter, 1^, which was defined as the fraction of Xi converted to 
C~ , is essentially the yield of biomass on Si. Indeed, substituting the 
relation, r^^i ~ Vs^i in (??) shows that after quasisteady state has been 
attained, approximates the fraction of Si that accumulates in the cell. 

(2) Although we have assumed that Yi is constant, the validity of (??) does 
not rest upon this assumption. It is true even if the yields are variable 
(i.e., depend on the physiological state). 

(3) Enzyme induction is subject to positive feedback. This becomes evident 
if (??) is substituted in (??): The quasisteady state induction rate, r^^, 
is an increasing function of e^. The positive feedback refiects the cyclic 
structure of enzyme synthesis that is inherent in Fig. ??: Ei promotes the 
synthesis of Xj, which, in turn, stimulates the induction of E^. 

We shall appeal to these facts later. 

The equations for batch growth are obtained from the above equations by let- 
ting D = 0. In what follows, we shall be particularly concerned with the initial 
evolution of the enzymes in experiments performed with sufficiently small in- 
ocula. Under this condition, the substrate consumption for the first few hours 
is so small that the substrate concentrations remain essentially constant at 
their initial levels, So,i. In the face of this quasiconstant environment, the en- 
zyme activities move toward a quasisteady state corresponding to exponential 
(balanced) growth. This motion is given by the equations 

de- 

= <i - iXg + Ki) ei, (13) 

Xi ~ , (^i4J 

i^x^i 

Tg ~ yiV;,iei , + Y^Vs,2e2 , , (lo) 

J^s,\ + ^0,1 J^s,2 -r So,2 

obtained by letting so,i in eqs. (??)-(??). The steady state(s) of these 
equations represent the quasisteady state enzyme activities attained during 
the first exponential growth phase. 

We note finally that the cells sense their environment through the ratio, 



K • + s- ' 



which characterizes the extent to which the permease is saturated with Si. 
Not surprisingly, we shall find later on that the physiological steady states 
depend on the ratio, CTj, rather than the substrate concentrations per se. It is, 
therefore, convenient to define the ratios 

= -p — — , = 
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For a given cell type, ai, ao,i. cT/^j are increasing functions of Sj, Si^, Sf^i, re- 
spectively, and can be treated as surrogates for the corresponding substrate 
concentrations. In what follows, we shall frequently refer to CTj, (To,i, and cr/^i 
as the substrate concentration, initial substrate concentration, and feed con- 
centration, respectively. 



3 Results 



3. 1 The role of regulation 

Before analyzing the foregoing equations for batch and continuous cultures, 
we pause to consider the role of regulatory mechanisms. This is important 
because the literature places great emphasis on the regulatory mechanisms, 
whereas the model neglects them completely (Vs,i and Ve^i are assumed to be 
constants). It is therefore relevant to ask: To what extent do the regulatory 
mechanisms control gene expression? To address this question, we shall focus 
on the well-characterized lac operon in E. coli. The analysis shows that cAMP 
activation and inducer exclusion cannot explain the strong repression observed 
in experiments. 



3.1.1 cAMP activation 

Interest in cAMP as a lac regulator began when it was discovered that the lac 
induction rate increased ~2-fold upon addition of cAMP to glucose-limited 
cultures (??). However, this remained a hypothesis until Epstein et al showed 
that the lac induction rate increased with the intracellular cAMP levels (Fig. ??a). 
They obtained this relationship by measuring the /3-galactosidase activities 
and intracellular cAMP levels during exponential growth of the /ac-constitutive 
strain E. coli X2950 on a wide variety of carbon sources. The use of a consti- 
tutive strain served to eliminate inducer exclusion: In such cells, the repressor- 
operator binding is impaired so severely that ctj ~ 0, and the induction rate 
is unaffected by the inducer level (r^^ ^ K,*)- The diverse carbon sources 
provided a way of varying the intracellular cAMP levels: In general, the lower 
the specific growth rate on a substrate, the higher the intracellular cAMP 
level (?). 

Although Fig. ??a shows a correlation between the cAMP level and the lac 
induction rate, it does not prove that there is a causal relationship between 
them. Indeed, since protein degradation is negligible in substrate-excess batch 
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Figure 6. The relative roles of dilution and cAMP activation in lac induction, 
(a) Variation of the /3-galactosidase activity with intracellular cAMP level during 
exponential growth of the /ac-constitutive strain E. coli X2950 on various carbon 
sources (?, Fig. 2). (b,c) During batch growth of (b) E. coli MC4100 ACPTIOO (?, 
Fig. 2) and (c) E. coli NC3 (?, Fig. 1) on various carbon sources, the steady state 
/3-galactosidase activity is inversely proportional to the specific growth rate on the 
carbon source. Moreover, the /3-galactosidase activity increases at most 2-fold if 
5mM cAMP is added to the medium (■). (d) The inverse relationship is satisfied 
at > 0.4 h"-*^ in continuous cultures of E. coli K12 W3110 grown on glucose and 
glucose + ImM IPTG (?, Fig. 1). The curves in (b)-(c) and (d) show the fits to 
(??) and (??), respectively. 



cultures (?), (??) yields 



r+. V ■ 

~:££~^. (16) 

^9 



Thus, the steady state enzyme level depends on the ratio of the specific induc- 
tion and growth rates. Since the substrates that yield high intracellular cAMP 
levels also support low specific growth rates, it is conceivable that Fig. ??a 
reflects the variation of the specific growth rate, Vg, rather than the specific 
induction rate, Ve^i- 

Eq. (??) provides a simple criterion for checking if the specific induction rate 
varied in the experiments: Evidently, Vg.j varies if and only if e.; is not inversely 
proportional to Vg. Unfortunately, we cannot check the data in Fig. ??a with 
this criterion, since the authors did not report the specific growth rates of the 
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carbon sources. However, the very same experiments have been performed by 
others (??), and this data shows that the /9-galactosidase activity is inversely 
proportional to Tg (Figs. ??b,c). It follows that in Figs. ??b and ??c, the 
galactosidase activity varies almost entirely due to dilution. The influence of 
cAMP on lac transcription is very modest, a fact that is further confirmed 
by the very marginal (< 2-fold) improvement of the /9-galactosidase activity in 
the presence of 5mM cAMP (Figs. ??c). The same is likely to be true of the 
data in Fig. ??a. Thus, the data from ?, which is extensively cited to support 
the existence of cAMP-mediated lac control, forces upon us just the opposite 
conclusion: cAMP has virtually no effect on the lac induction rate. 

The data from continuous cultures also implies the absence of significant 
cAMP activation. Indeed, (??) implies that when /ac-constitutive mutants 
of E. coli are grown in continuous cultures, the steady state /^-galactosidase 
activity at sufficiently high dilution rates is given by the relation 

e^ ^ (17) 

This relationship also holds for wild-type cells if the medium contains satu- 
rating concentrations of the gratuitous inducer, IPTG. In both cases, cAMP 
control is significant only if is not inversely proportional to D. But ex- 
periments performed with /ac-constitutive mutants (Fig. ??a) and wild-type 
cells exposed to saturating IPTG levels (Fig. ??d) show that /3-galactosidase 
activity is inversely proportional to D for sufficiently large D. ^ 

Catabolite repression also appears to be insignificant in cell types other than 
E. coli. The amidase activity is inversely proportional to D in the constitutive 
mutant P. aeroginosa Cll (Fig. ??a). Moreover, this declining trend is not due 
to catabolite repression, as postulated by Clarke and coworkers — it reflects 
the eff'ect of dilution. Likewise, the activity of alcohol oxidase is inversely 
proportional to D when H. polymorpha and C. boidinii are grown in glucose- 
limited chemostats. In this case, (??) implies that at sufficiently large dilution 
rates 

e.= T^^, (18) 
which is formally similar to (??). 



^ At small D, the enzyme activity is significantly lower than that predicted by (??). 
Evidently, this discrepancy cannot be due to cAMP activation. It cannot be resolved 
by enzyme degradation either: Although the data can be fitted to the equation, Cj = 
Ve^i/{D + ke,i), the values of fce,i required for these fits are physiologically implausible 
(~0.5 h^^). Based on the detailed study of the ptsG operon at various dilution rates, 
the discrepancy is probably due to enhanced expression of the starvation sigma 
factor, RpoS, at low dilution rates (?, Fig. 3) 
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3.1.2 Inducer exclusion 



The data shows that inducer exclusion certainly exists. When glucose or 
glucose-6-phosphate are added to a culture growing exponentially on lactose, 
the specific uptake rate of lactose declines ~2-fold (?, Figs. 3,5). Likewise, in 
the presence of glucose, the intracellular inducer concentration, as measured 
by the accumulation of radioactively labeled TMG, also decreases ~2-fold (?, 
Fig. 1). But this modest decrease cannot account for the dramatic repression 
of the lac operon during diauxic growth: The addition of glucose-6-phosphate 
to a culture growing on lactose decreases the activity of /?-galactosidase ~300- 
fold after only 4 generations of growth (?, Table 1). 

We conclude that cAMP activation and inducer exclusion do influence the lac 
induction rate, but their effect is not commensurate with the experimentally 
observed repression. We show below that the minimal model predicts complete 
transcriptional repression, despite the absence of any regulatory mechanism. 



3.2 Single-substrate growth 



In wild-type cells, the repressor-operator binding is so tight that ctj is large 
compared to 1. This is true even in the case of glucose (0:^=5-10), which is 
often treated as a substrate consumed by constitutive enzymes (?, Table 3). 
Under these conditions, the basal induction rate is small compared to the 
fully induced induction rate. Hence, at all but the smallest inducer levels, the 
specific induction rate can be approximated by the expression 

and the quasisteady state induction rate, obtained by substituting (??) in 
(??), is 

-"■6,1 ~r (^i^i ''s,i 

It is useful to make this approximation because the equations become amenable 
to rigorous analysis and yield simple physical insights. In the rest of the paper, 
we shall confine our attention to these idealized perfectly inducible systems. 

We begin by considering the growth of cultures limited by a single substrate, 
say. Si. Under this condition, S2 — €2 — 0, and the steady state values of the 
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Figure 7. Bifurcation diagram for single-substrate growth on Si. The blue curve 
denotes the curve of the critical dilution rate, -Dc,i, and denotes the threshold 
feed concentration. The 3-D insets depict the existence and stability of the steady 
states at any given ctj^i and D. Stable and unstable steady states are represented by 
full and open circles, respectively. The dashed line shows the transition dilution rate, 
Dt,i, which plays a crucial role in mixed-substrate growth (Section ??). The figure 
is not drawn to scale: Since /ce,i ~ 0.05 h~^ and Ve^i/Ke^i ~ 1 h~^, the threshold, 
al = A;e,i/(Ve,i/-^e,i), is on the order of 0.05. 



remaining variable satisfy the equations 



= ^ = - si) - K.ieiCTiC, (20) 

= ^ = Ve,i^r^ {YiVs,ieiai + ei, (21) 

at Ke,l + CiCTi 

dc 

= ^ = iY^Vs,lela^ - D) c, (22) 

X. - (23) 



These equations have three types of steady states: ei > 0, 62 = 0, c > 0, 
denoted -Eioi; ei > 0, 62 = 0, c = 0, denoted -Eioo; and ei = 0, 62 = 0, c = 0, 
denoted -Eqoo- The first two steady states correspond to the persistence and 
washout steady states of the classical Monod model (?). The third steady state, 
-E0007 is a consequence of perfect inducibility — it exists precisely because the 
specific induction rate is zero in the absence of the inducer. 



Fig. ?? shows the bifurcation diagram for eqs. (??)-(??). It can be inferred 
from the following facts derived in Appendix ??. 
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£^000 always exists, but it is stable if and only if 



Thus, growth cannot be sustained in a chemostat if the feed concentration 
is below the threshold level, a^. Under this condition 

Ke,l + ei(Ti Ke,l + ei(7/,i 

i.e., the induction rate of Ei never exceeds the degradation rate. Hence, 
as t ^ 00, ei approaches zero, and c,ai tend toward 0,a/,i, respectively. 
-Eioo exists if and only if ct/^i > a^, in which case it is given by the 
relations, c — 0, ai — cr/^i, and 



= ■ 

It is stable if and only if D exceeds the critical dilution rate 

Del = Y^Vs,ieiaf,i. (24) 

Evidently, £)c,i is an increasing function of (7/,i that is zero when cr/,i = al- 
We show below that L>c,i is the maximum specific growth rate that can 
be sustained in the chemostat (r^ < -Dc,i)- 

At first sight, this steady state poses a paradox: There are no cells in 
the chemostat, and yet, these non-existent cells have positive enzyme 
levels. The paradox disappears once it is recognized that £'100 reflects the 
long-run behavior of a small inoculum introduced into a sterile chemostat 
operating at ct/ i > u* and D > D^^i. As t — cxd, the specific growth rate 
of the cells approaches the maximum level, Dc^i. However, since Tg < £>c,i; 

dc 

— ^{rg-D)c< (L'ci - i^) c < 0, 

so that c — > 0, (Ti — > (T/,1, and ei approaches the positive value given 
above. Thus, in contrast to Eqqq, the cells fail to accumulate despite 
sustained induction of the enzymes because the dilution rate is higher 
than the maximum specific growth rate that can be attained. 
-E'loi exists if and only if ct/^i > al and D < Dc,i, in which case it is given 
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by the relations 



D 



(25) 



Ve,l D 



(26) 



L' + A;e,iL' + yiF.,iKe,i' 

D _ {d + k,,){d + y^v,,,k,,^) 



(27) 



c^Yi{sf,i- si), 



(28) 



all of which follow from the fact that the specific substrate uptake rate 
is proportional to D, i.e., 



Furthermore, i^ioi is stable whenever it exists. Evidently, this is the only 
steady state that allows positive cell densities to be sustained. 
The critical dilution rate, i?c,i, is the maximum specific growth rate that 
can be maintained at steady state. To see this, observe that at i?ioi, the 
specific growth rate equals the dilution rate. Thus, (??) says that specific 
growth rate increases with ai, and is maximal when cti = ctj^i. Letting 
(7i = in (??) and solving for D yields (??). 

Fig. (??) shows that at any given dilution rate and feed concentration, exactly 
one of the three steady states is stable. It is therefore straightforward to deduce 
the variation of the steady states along any path in the cx/^i, D-plane. The 
experiments are typically performed by varying either (T/,i or D, while the 
other one is held constant. We shall confine our attention to these two cases. 

3.2.1 Variation of steady states with at fixed D 

If the feed concentration is increased at a sufficiently small fixed D (horizontal 
arrow in Fig. ??), steady growth is feasible only if the feed concentration lies 
in the region to the right of the blue curve, where i^ioi is stable. The variation 
of the steady states in this region is given by (??)-(??), which imply that 
as the feed concentration increases, xi, ei, cti remain constant, and c increases 
linearly. Thus, any increase in the feed concentration is compensated by a pro- 
portional increase in the cell density — all other variables remain unchanged. 
We shall appeal to this fact later. 

No data is available for the inducer or enzyme levels, but experiments per- 
formed with glucose-limited cultures of E. coli ML308 show that at dilution 
rates of 0.3 and 0.6 h^\ the residual substrate concentration is indeed inde- 
pendent of the feed concentration, and the cell density does increase linearly 
with the feed concentration (Fig. ??a). No threshold was observed in these ex- 



D 



(29) 
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Figure 8. Growth of E. coli ML308 in a glucose-limited chemostat. (a) When the 
feed concentration of glucose is increased at a fixed dilution rate, the residual glu- 
cose concentration (closed symbols) is constant, and the cell density (open symbols) 
increases linearly (?, Appendix B). The data for D = 0.3 h~^ and D = 0.6 h~^ 
is represented by the symbols. A, A and respectively, (a) When the dilution 

rate is increased at a fixed feed concentration (100 mg L~^), the residual substrate 
concentration increases monotonically with D, and approaches a threshold at small 
dilution rates (?, Fig. 2). The data was obtained at 4 temperatures. 

periments because the feed concentrations were relatively high (> 1 mg L~^). 
As we show below, the threshold concentration for glucose is ~10-50 /ig L~^. 



3.2.2 Variation of steady states with D at fixed (7/,i 

If the dilution rate is increased at any fixed cr/^i < a^, there is no growth 
regardless of the dilution rate at which the chemostat is operated. If uj^i > a^, 
the variation of the steady states with D is qualitatively identical to that of the 
Monod model: The persistence steady state, £'ioi, is stable for < D < i?c,i, 
and the washout steady state, £"1005 is stable for D > Dc,i (vertical arrow in 
Fig. ??). We shall confine our attention to £^101, since £100 is independent of 
D. 



Variation of inducer and enzyme levels It follows from (??) (??) that 
as the dilution rate increases over the range, (0,Dc,i), the inducer level in- 
creases linearly, whereas the enzyme level passes through a maximum at 

D — y ^YiV^^i^e,i) ke,i- The latter is consistent with the data for systems 
with inducible enzymes (Fig. ??). 

In contrast to the hypothesis of Clarke et al., the model explains the non- 
monotonic profile of the enzyme level as the outcome of a balance between 
induction and dilution/degradation (as opposed to catabolite repression). To 
see this, observe that the steady state enzyme level is given by the ratio of the 
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induction and dilution/degradation rates, i.e., 



Furthermore, since xi increases with D, so does r^^, but it saturates at suffi- 
ciently large D. At low values of the enzyme level increases with D because 
the induction rate increases with D. At high D, the steady state enzyme level 
inversely proportional to D because the induction rate saturates. 

The validity of the above explanation rests upon the hypothesis that the in- 
duction rate saturates at large D. This hypothesis is supported by the data 
for the lac and ptsG operons. When batch cultures of E. coli growing ex- 
ponentially on saturating concentrations of lactose (which are physiologically 
identical to cultures growing in a lactose-limited chemostat near the critical 
dilution rate) are exposed to 0.4 mM IPTG, there is almost no change in the 
/3-galactosidase activity (?, Table 1). Similarly, the PtsG activity of wild-type 
and pisG-constitutive E. coli is the same in glucose-excess batch cultures and 
continuous cultures operating at high dilution rates (?, Fig. 3). Induction is 
therefore already saturated during growth of E. coli on lactose or glucose, 
provided the dilution rate is sufficiently large. It remains to be seen if this 
hypothesis is also valid for the other systems. However, at large dilution rates, 
the enzyme activity is more or less inversely proportional to D for all the sys- 
tems shown Fig. ??, a result expected only if the induction rate saturates at 
large D. Thus, the decline of the enzyme activity in both constitutive mutants 
and wild-type cells is due to dilution, rather than catabolite repression. 



Variation of the substrate concentration and cell density According 
to the Monod model, the residual substrate concentration is proportional to 
the dilution rate, i.e., ai = D /r™f^, where r™f^ denotes the maximum specific 
growth rate on Si (?). In contrast, eq. (??) yields 

which implies that at small dilution rates, iTi approaches the threshold level, 
al, and at large dilution rates, cti is proportional to D^. Both these properties 
follow from the relation, r^^ = Vg^ieiai = D/Yi, and the steady state profiles 
of the peripheral enzyme activity. Indeed, at low dilution rates, ai is constant 
because the enzyme level increases linearly. Likewise, at high dilution rates, 
(Ti ~ D'^ because ei ~ 1/D. 

The experiments provide clear evidence of threshold concentrations. Kovarova 
et al. have shown that when E. coli ML308 is grown in a glucose-limited 
chemostat, the residual glucose concentration approaches threshold levels of 
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20-40 /xg/L at small dilution rates (Fig. ??b). The authors did not identify 
the specific mechanism imderlying the threshold, ft may exist, for example, 
simply because the maintenance requirements preclude growth. However, we 
shall show below that in some cases, the threshold exists precisely because the 
enzymes are not induced at sufficiently low values of the substrate concentra- 
tion. 



3.3 Mixed- substrate growth 



We have shown above that in chemostats limited by a single substrate, enzyme 
synthesis is abolished if the feed concentration is below the threshold level, 
and this occurs because the induction rate of the enzyme cannot match the 
degradation rate. Our next goal is to show that in mixed-substrate cultures, 
enzyme synthesis is often abolished because the induction rate of the enzyme 
cannot match its dilution rate. To this end, we shall begin with the case of 
mixed-substrate batch cultures, where this phenomenon is particularly trans- 
parent. This analysis of batch cultures will also enable us to establish the 
formal similarity between the classification of the growth patterns in batch 
and continuous cultures. 



3.3.1 Batch cultures 

3.3.1.1 Mathematical classification of the growth patterns The 

growth of mixed-substrate batch cultures is empirically classified based on 
the substrate consumption pattern (preferential or simultaneous) during the 
first exponential growth phase. The mathematical correlate of the foregoing 
empirical classification is provided by the bifurcation diagram for the equations 

^ = Ve,i^^^ {YiVs,ieiao,i + y2K,2e2<7o,2 + ke,i) ei, (30) 

dt Ag^i + 61(70,1 

= Ve,2^r-^^^ (i^lK,ieiCro,i + 5^2^^5,2620-0,2 + ke,2) 62, (31) 

dt Ae,2 + 62(70,2 

which describe the motion of the enzymes toward the quasisteady state cor- 
responding to the first exponential growth phase. 

We gain physical insight into the foregoing equations by observing that they 
are formally similar to the Lotka-Volterra model for two competing species 

^ = a^N, - {bM + b^2N2) N,, (32) 
dt 

^ = a2N2 - {b2iN, + b22N2) N2, (33) 
dt 
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where Ni is the population density of the i-th species; OjA^i is the (unrestricted) 
growth rate of the ?'-th species in the absence of any competion: and buNf ^ 
bijNiNj represent the growth rate reduction due to intra- and inter-specific 
competition, respectively. Evidently, the net induction rate, Ve^ieiao^i/{Ke^i + 
^iCo,i) — kciSi, and the two dilution terms in (??)-(??) are the correlates 
of the unrestricted growth rate and the two competition terms in the Lotka- 
Volterra model. Thus, the interaction between the enzymes is competitive (the 
enzymes inhibit each other via the dilution terms), despite the absence of any 
transcriptional repression (the induction rate of Ei is unaffected by activity of 

Now, the Lotka-Volterra equations are known to yield coexistence or extinc- 
tion steady states, depending on the values of the parameters, bij (?). The 
formal similarity of eqs. (??)-(??) and (??)-(??) suggests that the enzymes 
will exhibit coexistence and extinction steady states, depending on the values 
of the physiological parameters and the initial substrate concentrations. We 
show below that this is indeed the case. Moreover, the coexistence and extinc- 
tion steady states are the correlates of simultaneous and preferential growth, 
respectively. 

Eqs. (??)-(??) have 4 types of steady states: ei = 0, 62 = 0; ei > 0, 62 = 0; 
ei = 0, 62 > 0; and ei > 0, 62 > 0, which are denoted Eqq, Eiq, Eqi, and En, 
respectively. Each of these steady states has a simple biological interpretation. 
If the enzyme levels reach £'00 (resp., Eu), none (resp., both) of the two 
substrates are consumed during the first exponential growth phase. If the 
enzyme levels reach £"10 (resp., £01), only Si (resp., S2) is consumed during the 
first exponential growth phase. The bifurcation diagram shows the existence 
and stability of these steady states (and hence, the growth pattern) at any 
given (To,i and (7o,2 (Fig- ??)■ It can be inferred from the following facts derived 
in Appendix ??. 

(1) £00 always exists (for all (To,i and (70^2) ■ It is stable if and only if 

r;,i = - ke, < ^ ao, < a*, t = 1,2 

i.e., the point, (cro,i, cro,2), lies below the green line and to the left of the 
black line in Fig. ??. 

(2) £10 exists if and only if (Jqi > a^, in which case it is unique, and given 

by 



-(^e, + ^) + V(^e, + ^) + 



ei = , 62 = 0. 

At this steady state, the cells consume only 5*1, and grow exponentially 
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Figure 9. Bifurcation diagram classifying the growth patterns of mixed-substrate 
batch cultures at various initial substrate concentrations, cro,i and cjo,2- The green 
and black lines represent the threshold concentrations for induction of Ei and E2 , 
respectively. The blue and red curves represent the locus of initial substrate con- 
centrations defined by the equations, r^^i = r*2 and rg^2 = respectively. The 
insets in each of the six regions show the dynamics at the corresponding initial sub- 
strate concentrations. The first inset shows the motion of the enzymes during the 
first exponential growth phase; the second one shows the resultant evolution of the 
substrate concentrations and cell density. 

at the specific growth rate 



rg,i = FiFs.ieicro,!. 



It is stable if and only if 

< rg,i, (34) 
i.e., (cTo,i, cro,2) lies below the blue curve in Fig. ??, defined by the equation 

= rg,i. (35) 

Thus, £'10 is stable under two distinct sets of conditions. Either (T2,o is 
below the threshold level, cTg, in which case synthesis of E2 cannot be 
sustained because its induction rate cannot match its degradation rate. 
Alternatively, (To,2 is above the threshold level, but sustained synthesis 
of E2 remains infeasible because the dilution of E2 due to growth on Si 
is too large. Thus, r* 2 can be viewed as the highest specific growth rate 
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tolerated by the induction machinery for E2: Synthesis of E2 is abolished 
whenever r^^i exceeds r*2- We shall refer to r*2 as the specific growth 
rate for extinction of E2. 

For any given ctq,! > cr^, the exponential growth rate, r^^i, can be higher or 
lower than r* 2 (depending on the value of (70,2), but it cannot exceed the 
corresponding extinction growth rate, r* i. This reflects the biologically 
plausible fact that ei cannot be zero during exponential growth on 5*1. 
(3) Eqi exists if and only if (To,2 > (^2^ which case it is (uniquely) given by 
the expressions 



ei = 0, 62 = 



2 4Ke 2r 



2cro,2 

At this steady state, the cells consume only 5*2, and grow exponentially 
at the specific growth rate 

rg,2 = ^2^8,2620-0,2 ■ 

It is stable if and only if 

r*g,i < rg,2, (36) 
i.e., ((7o,i, (70,2) lies to the left of the red curve in Fig. ??, defined by the 
equation 

rli = rg,2. (37) 
Furthermore, 2 < r* 2 with equality being attained precisely when (72 = 
(72, in which case rg^2 = ''^^,2 = 0- 

(4) £'11 exists and is unique if and only if £'10 and £01 are unstable, a condition 
satisfied precisely when 

i.e., (ctq^i, (70,2) lie between the blue and red curves in Fig. ??. At this 
steady state, both enzymes "coexist" because neither substrate supports 
a specific growth rate that is large enough to annihilate the enzymes for 
the other substrate. Hence, both substrates are consumed, and the cells 
grow exponentially at the specific growth rate 

rg,i2 = i^iK,iei<7o,i + ^2^5,262 0-0,2, 

where 61,62 are the unique positive solutions of (??)-(??). 

(5) The blue curve and red curves are given by the equations 



K..2 \ . 1 



Y2Vs,2K,,2j Y2Vs,2K,,2 
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respectively. Both equations define graphs of increasing functions passing 
through (ci,C2), but the blue curve always lies below the red curve. 

Thus, the bifurcation diagram consists of 6 distinct regions such that exactly 
one of the four steady states is stable in each region. 

The bifurcation diagram implies that there are six possible growth patterns, 
depending on the initial substrate concentrations. Three of these growth pat- 
terns occur if (To,i < (7* or/and (To,i < (72. If both a"o,i and (7o,2 are below their 
respective threshold concentrations, neither substrate is consumed. If only one 
of them is below its threshold level, say, (To,2 < o'2 5 '^^h' is consumed during 
the first exponential growth phase. However, growth is not diauxic because 5*2 
is never consumed. If ao,i > crl and (To,2 > C2, the model predicts the existence 
of three additional growth patterns. 

(1) If (To,i, (70,2 lie in the region between the blue and green curves, 62 ap- 
proaches zero during the first exponential growth phase. Importantly, 
since (To,2 is higher than the threshold level, E2 is synthesized upon ex- 
haustion of 5*1 at the end of the first exponential growth phase. Hence, 
growth is diauxic with preferential consumption of Si. 

(2) If (To,i, (^0,2 lie in the region between the red and black curves, ei ap- 
proaches zero during the first exponential growth phase, i.e., there is 
diauxic growth with preferential consumption of S2. 

(3) If ctq,!, (7o,2 lie between the red and blue curves, Ci and 62 attain posi- 
tive values during the first exponential growth phase. Consequently, both 
substrates are consumed until one of them is exhausted, and there is no 
diauxic lag before the remaining substrate is consumed. 

We show below that the foregoing classification is consistent with the data. 



3.3.1.2 Growth patterns at saturating substrate concentrations 

Physiological experiments, which are generally performed with saturating sub- 
strate concentrations ((7o,i,(7o,2 ~ 1), show that depending on the cell type, 
the two substrates are consumed sequentially or simultaneously. The behav- 
ior of the model is consistent with this observation. To see this, observe that 
depending on the values of the physiological parameters (which, in effect, de- 
termine the cell type), there are three possible arrangements of the red and 
blue curves, each of which yields a distinct growth pattern at saturating sub- 
strate concentrations. Indeed, Si is consumed preferentially if (1, 1) lies below 
the blue curve (Fig. ??a). This occurs precisely when the physiological param- 
eters satisfy the condition 




(40) 
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Figure 10. Bifurcation diagrams corresponding to diauxic growth under substrate-ex- 
cess conditions (cTo,i,cro,2 ~ !)• (a) Preferential consumption of Si. (b) Preferential 
consumption of 82- 

obtained from (??) by neglecting enzyme degradation (ke,i = 0). Likewise, 5*2 
is consumed preferentially if (1, 1) lies above the red curve (Fig. ??b), i.e., the 
physiological parameters satisfy the condition 

^-->Sl-,^^^/^^ fSlV, (41) 



Ke,2 Ke,l Y2Vs,2Ke,U \Ke,l, 

obtained from (??) upon setting fce,i = 0. Finally, the substrates are consumed 
simultaneously if (1,1) lies above the blue curve and below the red curve 
(Fig. ??), i.e., both the foregoing conditions are violated. 

The biological meaning of these conditions was discussed in earlier work (?). 
Indeed, one can check that (??) and (??) are equivalent to the conditions 



a < 
a > a* 



2 
2 

(^-K2 + y^Kf + 4^ 



respectively, where a = ^^12^,2^,2/ y^^iK.iK,! is a measure of the maximum 

specific growth rate on S2 relative to that on 5*1, and Ki = K^^i/ ^jVeJJiXiVsJ) 
is a measure of saturation constant for induction of Ei. Thus, the substrates 
are consumed preferentially whenever the maximum specific growth rate on 
one of the substrates is sufficiently large compared to that on the other sub- 
strate. Just how large this ratio must be depends on the saturation constants, 
Ki. Enzymes with small saturation constants are quasi-constitutive — their 
synthesis cannot be abolished even if the other substrate supports a large 
specific growth rate. In contrast, synthesis of enzymes with large saturation 
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Figure 11. Transitions of the growth patterns triggered by changes in the initial 
substrate concentrations. Upper panel: Batch growth of Escherichia coli ML308 
on mixtures of glucose and galactose, (a) At high initial glucose concentrations 
(20 mg L-i), glue ;ose is consumed before galactose (?). (b) At low initial glucose 
concentrations (5 mg L^^), glucose and galactose are consumed simultaneously (?, 
Fig. 1). Lower panel: Batch growth of Pseudomonas on PNP and PNP + glu- 
cose (?, Figs. 1,4). DPM denotes disintegrations per min of -"^^C-labelled PNP and 
glucose, (c) No PNP is consumed if its initial concentration is 10 ^g L~^. Significant 
consumption occurs at initial concentrations of 50 and 100 //g L~^. (d) PNP and 
glucose are consumed simultaneously if their initial concentrations are high (3 and 
10 mg L~^, respectively). 

constants can be abolished even if the other substrate supports a comparable 
specific growth rate. 



3.3.1.3 Transitions triggered by changes in substrate concentra- 
tions Figs. ??-?? imply that the growth pattern can be changed by altering 
the initial substrate concentrations. We show below that this conclusion is 
consistent with the data. 

As a first example, consider a cell type that prefers to consume Si under 
substrate-excess conditions ((To,i,(To,2 ~ !)• Its growth pattern at any (To,i, 
(70,2 is then described by Fig. ??a, w^hich implies that if ao.i is decreased 
sufficiently from 1, both substrates will be consumed simultaneously. Such 
behavior has been observed in experiments. If the initial concentrations of 
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glucose and galactose are 20 mg/L and 5 mg/L, respectively, E. coli ML308 
consumes glucose before galactose (Fig. ??a). If the initial concentration of 
glucose is reduced to 5 mg/L or lower, glucose and galactose are consumed 
simultaneously (Fig. ??b). 

The data in Figs. ??a,b are open to question because the initial cell densities 
are so large (~ 10^ cells L^^) that the substrates may have been exhausted 
before the enzymes reached the quasisteady state corresponding to balanced 
growth. However, similar results have also been obtained in studies performed 
with very small initial cell densities (10^-10^ cells L~^). Indeed, Schmidt and 
Alexander observed that "simultaneous use of two substrates is concentration 
dependent . . . Pseudomonas sp. ANL 50 mineralized glucose and aniline si- 
multaneously when present at 3 fig but metabolized them diauxically at 
300 fig L~^" (?, Figs. 3-4). Likewise, Pseudomonas acidovorans consumed ac- 
etate before phenol when the concentrations of these substrates were >70 and 
2 yug L~^, respectively. If the initial concentration of acetate was reduced to 
13 /ig L~^, both substrates were consumed simultaneously. 

In terms of the model, these transitions from sequential to simultaneous sub- 
strate consumption can be understood as follows. Sequential substrate con- 
sumption occurs at high concentrations of the preferred substrate because it 
supports such a large specific growth rate that the enzymes associated with 
the secondary substrate are diluted to extinction. Reducing the initial con- 
centration of the preferred substrate decreases its ability to support growth, 
and thus dilute the secondary substrate enzymes to extinction. Consequently, 
both substrates are consumed. 

The growth of Pseudomonas sp. on p-nitrophenol (PNP) and glucose furnishes 
additional examples of growth pattern transitions driven by changes in the ini- 
tial substrate concentrations. Schmidt et al. found that no PNP was consumed 
if its initial concentration was 10 /ig L~^, and significant mineralization oc- 
curred at the higher initial concentrations of 50 and 100 fig (Fig. ??c). 
It follows that the threshold concentration for PNP lies between 10 and 50 
fig L~^. Importantly, the addition of 20 mg L~^ glucose to a culture containing 
10 fj,g PNP failed to stimulate PNP consumption. The authors concluded 
that "the threshold was not a result of the fact that the concentration of PNP 
was too low to meet maintenance energy requirements of the organism, but 
rather supports the hypothesis that the concentration was too low to induce 
degradative enzymes." Finally, it was observed that when the initial concen- 
trations of PNP and glucose were increased to saturating levels of 3 mg L~^ 
and 10 mg respectively, both substrates were consumed simultaneously 
(Fig. ??d). All these results are consistent with the bifurcation diagram shown 
in Fig. ??. 

Figs. 11 11 extend the results obtained in ?. There, we constructed bifurca- 
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tion diagrams describing the variation of the growth patterns in response to 
changes in the physiological parameters at (fixed) saturating substrate con- 
centrations. These genetic bifurcation diagrams explained the phenotypes of 
many different mutants and recombinant cells. Figs. ??-?? can be viewed 
as epigenetic bifurcation diagrams, since they describe the variation of the 
growth pattern when a given cell type, characterized by fixed physiological 
parameters, is exposed to various initial substrate concentrations. 



3.3.2 Continuous cultures 

In the experimental literature, the growth pattern of mixed-substrate contin- 
uous cultures is classified based on the manner in which the two substrates 
are consumed when the dilution rate is increased at a given pair of feed con- 
centrations. Either both substrates are consumed at all dilution rates up to 
washout, or both substrates are consumed up to an intermediate dilution rate 
beyond which only one of the substrates is consumed. We show below that the 
model predicts these growth patterns. Moreover, it can quantitatively capture 
the variation of the steady state enzyme levels, substrate concentrations, and 
cell density with respect to the dilution rate and the feed concentrations. 



3.3.2.1 Mathematical classification of the growth patterns The 

mathematical correlate of the experimentally observed growth patterns cor- 
responds to the bifurcation diagram for eqs. (??)-(??), which describes the 
existence and stability of the steady states in the 3-dimensional cr/^i, -D- 
space. 

Eqs. (??)-(??) have 7 steady states, but 5 of them are "single-substrate steady 
states" (see rows 1-5 of Table ??). Indeed, Eqqq, -Eioi, and Eiqq (resp., i?ooo, 
i^oii, and Eqio) are observed in chemostats limited by 5*1 (resp., 5*2). Their 
existence in mixed-substrate growth implies that even if both substrates are 
supplied to the chemostat, there are steady states at which only one of the 
substrates is consumed. We shall see below that these steady states can become 
stable under certain conditions. Here, it suffices to observe that only two of 7 
steady states, namely, £"111 (ei > 0, 62 > 0, c > 0) and Ehq (ei > 0, 62 > 0, 
c > 0), are uniquely associated with mixed-substrate growth. They correspond 
to simultaneous consumption of both substrates and its washout. 

The existence and stability of the steady states are completely determined 
by five special dilution rates, namely, Dt^i, Dc,i, and Dc (Table 1, columns 4- 
5). These dilution rates are the analogs of the special specific growth rates 
considered in the analysis of batch cultures. Indeed, the transition dilution 
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Table 1 

Necessary and sufficient conditions for existence and stability of the steady states 
for mixed-substrate growth in continuous cultures (see Appendix C for details). 
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(42) 



is the analog of the extinction growth rate, r* j. It is the maximum dilution 
rate up to which synthesis of can be sustained. The critical dilution rate, 
A,i; given by the relation 



-(^e,+ ^) + V(i^e. + ^) + 



is the analog of rg^i. It is the maximum dilution rate up to which growth can 
be sustained in a chemostat supplied with only Si, and satisfies the relation, 
Dei < Dt^i with equality being obtained precisely when (Xj = a* , in which case, 
A,j = A,i = (see dashed and blue lines in Fig. ??). Finally, the critical 



dilution rate. A, is defined as 



A = yiVs,ieiaf^i + Y2Vs,2e2crf,2: 
where 61,62 are the unique positive steady state solutions of (??)-(??) with 
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Figure 12. Bifurcation diagram classifying the growth patterns of mixed-substrate 
continuous cultures at various feed concentrations, aj^i and The green and black 
lines represent the threshold concentrations for induction of Ei and E2, respectively. 
The blue and red curves represent the locus of initial substrate concentrations defined 
by the equations, L'c,! = Dt^2 and -Dc,2 = Df i, respectively. The insets in each of 
the six regions show the variation of the steady states with D at the corresponding 
feed concentrations. 



cro,i replaced by a/^j. It is the analog of the function, Vg^^. Thus, r* j, j, Tg^u, 
and their analogs, Dt^i, Dc^i, Dc, are defined by formally identical expressions, 
the only difference being that (To,i is replaced by a/^j. 

The formal similarity of the expressions for Dt i, D^^i^ and r* j, j, 12 
implies that for a given set of physiological parameters, the relations, u/^j = 
0"*, Dt^i = Dc,2 and Dt^ = -Dc,i, define curves on the cr/^i, (j/,2-plane that 
are identical to the curves on the (7o,i, o"o,2-plane defined by the equations, 
o'o.i = ^3,1 = rg^2 and r* 2 = r^^i. It follows that the relations, cr/,i = cr*, 
Df i = Dc^2, Dt 2 = -Dc,i, determine a partitioning of the o"/ 1, (Tj2-plane into 
6 regions, which is identical to the partitioning of the bifurcation diagram for 
batch cultures, the only difference being that the coordinates are crji,(T/2, 
rather than cro,i,cro,2 (Fig- 
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Each of the six regions in Fig. ?? corresponds to unique growth pattern that 
can be inferred from the existence and stability properties summarized in 
columns 4-5 of Table 1. For instance, if (cr/,i, cr/_2) lies in the region between 
the green and blue curves in Fig. ??, both substrates are consumed for all 
< D < Dt^2, only is consumed for all Dt^2 < -D < Dc,i, and neither 
substrate is consumed because the cells are washed out for all D > -Dc,i- To 
see this, observe that in the region of interest (between the green and blue 
curves), the special dilution rates stand in the relation 

< L'e,2 < A,2 <D,< L>e,l < A,l- 

The existence and stability conditions listed in Table ?? therefore imply that 
as D increases, 4 of the 7 steady states play no role because they are unstable 
whenever they exist, or do not exist at all. Indeed: 

(1) Eqqq and Eqiq exist at all D > 0, but they are always unstable because 
the stability conditions for these two steady states {Dt^i,Dt2 > and 
Df i < Dc^2, respectively) are violated in the region of interest . 

(2) Eon exists for all < D < -Dc,2, but it is unstable whenever it exists 
because the stability condition, D > i, is not satisfied in the region of 
interest. 

(3) -Eiio does not exist at all since one of the existence conditions, Dt 2 > -Dc,i, 
is violated in the region of interest. 

It follows that only the remaining three steady states are relevant. One can 
check by appealing to Table 1 that is stable for < < Dt^2, -^'loi is 
stable for Dt 2 < D < L>c,i5 and £"100 is stable for D > Dc,i- Hence, both 
substrates are consumed for all < D < Dt^2- At D — Dt 2, 62 becomes zero. 
Only 5*1 is consumed for all 2 < D < Dc,i until c becomes zero at D = -Dc,i- 
Neither substrate is consumed for D > Dc,i. 

Similar arguments, applied to all the regions in Fig. ??, show that there are 6 
distinct growth patterns. Three of these growth patterns occur when at least 
one of the feed concentrations is below its threshold level. 

(1) If both feed concentrations are below the threshold level, £^000 is always 
stable, so that neither substrate is consumed at any dilution rate. 

(2) If (7f^2 < o'25 but (Tf i > al, EiQi is stable (i.e., only is consumed) for 
all < D < -Dc,i5 and £100 is stable (i.e., the cells wash out) whenever 
D > -Dc,i- This is qualitatively similar to the behavior observed during 
single-substrate growth on Si (Fig. ??). Thus, if the feed concentration 
of S2 is at sub-threshold levels, it has no effect on the behavior of the 
chemostat. 

(3) If af i < al, and (7/2 > o"!' -^011 is stable (i.e., only 5*2 is consumed) 
for all < D < -Dc,2, and £010 is stable (i.e., the cells wash out) for 
D > Dc,2- This is qualitatively similar to the behavior of the chemostat 
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during single-substrate growth on S2. 

If both feed concentrations are above their respective threshold levels, there 
are three additional growth patterns. 

(1) If (7/,i and (j/,2 lie between the red and blue curves, Em (i.e., both sub- 
strates are consumed) foi < D < Dc, and Euo is stable (i.e., cells wash 
out) for D > Df,. Thus, both substrates are consumed at all dilution rates 
up to washout. 

(2) If (T/,1 and (T/,2 lie between the green and blue curves, Em is stable (i.e., 
both substrates are consumed) for a\\ < D < Dt2, -Eioi is stable (i.e., 
only Si is consumed) for all Dt^i < D < Dc,i, and £^100 is stable (the cells 
wash out) for all D > -Dc,i- 

(3) If CT/j and (J/_2 lie between the black and red curves, Em is stable (i.e., 
both substrates are consumed) for all < < 1, Equ is stable (i.e., 
only ^2 is consumed) for all 2 < D < L'c,2, and ii^ioo is stable (the cells 
wash out) for all D > £)c,2- 

These growth patterns are the mathematical correlates of the growth patterns 
shown in Fig. ??. 

The model predicts that the growth pattern in continuous cultures can be 
changed by altering the feed concentrations. Unfortunately, there is no experi- 
mental data at subsaturating feed concentrations (cr/^j < ^s,?) due to technical 
difficulties associated with wall growth. Henceforth, we shall confine our at- 
tention to growth at saturating feed concentrations (cr/^j ~ 1). 

We have shown above that for a given cell type, the bifurcation diagrams for 
batch and continuous cultures are formally identical — one can be generated 
from the other by a mere relabeling of the coordinate axes. This formal iden- 
tity provides a precise mathematical explanation for the empirically observed 
correlation between the batch and continuous growth patterns of a given cell 
type. To see this, suppose the cells in question consume preferentially in 
substrate-excess batch cultures (cro,i ~ 1). The bifurcation diagram for this 
system then has the form shown in Fig. ??a. The continuous growth patterns 
of this system are given by the very same figure, the only difference being that 
the coordinates, (To,i, (To,2, are replaced by cji, af^2, respectively. It follows 
that if the cells are grown in continuous cultures fed with saturating substrate 
concentrations (cr/^j ~ 1), they will consume both substrates at low dilution 
rates, and only Si at all high dilution rates up to washout (Fig. ??b). A similar 
argument shows that if the cells consume both substrates in substrate-excess 
batch culture (bifurcation diagram given by a relabeled Fig. ??), their growth 
in continuous cultures will be such that both substrates are consumed at all 
dilution rates up to washout (Fig. ??a). 
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3.3.2.2 The model predicts the observed variation of the steady 
states with D and Sf_i Although the model predicts the existence of the 
observed growth patterns, it remains to determine if the variation of the steady 
states with D and is in quantitative agreement with the data. Since the 
variation of the "single-substrate steady states" with D and Sf^i is consistent 
with the data (Section ??), it suffices to focus on £"111, the only growth- 
supporting steady state uniquely associated with mixed-substrate growth. 
This steady state can be observed in systems exhibiting both the simultaneous 
and the preferential growth patterns. However, since the data for simultaneous 
systems is limited, we shall focus on preferential systems. In our discussion of 
these systems, we shall assume, furthermore, that 5*1 is the "preferred" sub- 
strate. This entails no loss of generality since the equations for and S2 (and 
the associated physiological variables) are formally identical: Interchanging 
the indices does not change the form of the equations. 

Two types of experiments can be found in the literature. 

(1) Either the dilution rate was changed at fixed feed concentrations. 

(2) Alternatively, the mass fraction of the substrates in the feed, V^i = Sf^/ (s/_i-|- 

was changed at fixed total feed concentration, Sf^t = Sf,i + s/,2, and 
(sufficiently small) dilution rate. 

The model predicts that in the first case, Em is stable at all dilution rates 
below Dt,2, the transition dilution rate for the "less preferred" substrate. In the 
second case, Em is stable at all feed fractions, except those near the extreme 
values, and 1, where the model predicts threshold effects. Such threshold 
effects have been observed when the feed contained a very small fraction of 
one of the substrates (???). Nevertheless, we shall restrict our attention to 
intermediate values of ipi at which is the stable steady state. 

The steady state, £"111, satisfies the equations 

= D{sf^i - Si) - Ts^iC, (43) 

= Ve,J^ {D + ke,) , (44) 

= Fir,,i + Y2rs,2 - D, (45) 
Xi ^ (46) 

These equations do not yield an analytical solution valid for all D and s/^j. 
However, we show below that they can be solved for the two limiting cases of 
low and high D, from which the entire solution can be pieced together. 

Before considering these limiting cases, it is useful to note that eqs. (??) and 
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(??) imply that 



c = Yi - si) + Y2 (s/,2 - S2) , (47) 

D f"'- Y^isf,^-s,) + Y2isf,2-S2y ^ ' 

The first relation says that the steady state cell density is the sum of the 
cell densities derived from the two substrates. The second relation states that 
f3i, the fraction of biomass derived from Si, equals the fraction of the total 
specific growth rate, D, supported by S^. Evidently, A < 1 with equality 
being attained only during single-substrate growth. 

Eq. (??) immediately yields 

rs,-^, (49) 

which explains the following well-known empirical observation: At any given 
D, the mixed-substrate specific uptake rate of Si is always lower than the 
single-substrate specific uptake rate of 5'j (reviewed in ?). This correlation is a 
natural consequence of the mass balances for the substrates and cells, together 
with the constant yield property. The mixed-substrate specific uptake rate is 
smaller simply because both substrates are consumed to support the specific 
growth rate, D, whereas only one substrate supports the very same specific 
growth rate during single-substrate growth. 



3.3.2.3 Low dilution rates The first limiting case corresponds to di- 
lution rates so small that both substrates are at subsaturating levels, i.e., 
~ ^ Since both substrates are almost completely consumed, 
(??)-(??) become 

c^YiSf^i + Y2Sf,2, (50) 

piD YiSf^i Yitpi 
r,,- = , 3j a; = . (51) 

' rr ^ YiSf,i+Y2Sf,2 yiV'l+1^2V'2 ^ ^ 

Evidently, /3j is completely determined by ipi, the fraction of Si in the feed. 
Moreover, Pi{ipi) is an increasing function of ■^j with /3i(0) = and /3i(l) = 1. 
It is identical to ipi when the yields on both substrates are the same. 

Eq. (??) implies that if D is increased at fixed feed concentrations, rs,i increases 
linearly. Likewise, if ipi, and hence, is increased at a fixed dilution rate, 
rs,i increases linearly with /^j. Both conclusions are consistent with the data. 
Fig. ?? shows that eq. (??) provides excellent fits to the data for growth of 
H. polyrnorpha on glucose + methanol and E. coli on glucose I lactose. The 
validity of (??) at low dilution rates has been demonstrated for many different 
systems (reviewed in ?). 
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Figure 13. Upper panel: The specific substrate uptake rate increases linearly when 
the dilution rate is increased at fixed feed concentrations. Specific uptake rates of 
(a) glucose and (b) methanol during growth of H. polymorpha on various mixtures of 
glucose and methanol (?, Figs. 4-5). The percentages show the percent methanol in 
the feed. The total feed concentration of glucose and methanol was 5 g L~^. Lower 
panel: If the fraction of Si in the feed is increased at a fixed dilution rate, the specific 
substrate uptake rate increases linearly with Pi. (c) Growth of H. polymorpha on a 
mixture of glucose and methanol at D = 0.145 h~^ (?, Fig. 5). (d) Growth of E. 
coli K12 on various mixtures of lactose and glucose at D = 0.4 h~^ (calculated from 
Fig. 4 of ?). The lines in the figures show the specific substrate uptake rates predicted 
by eq. (??) with the experimentally measured single-substrate yields, Yglu = 0.55, 
y\iET = 0.38 in (a)-(c), and Fglu = >lac in (d). 

The physiological variables, Xi and Cj, are completely determined by the spe- 
cific uptake rate of Si. Indeed, (??) and (??) imply that 

" = ^^^^ 

' D+ke,iYiVs,iKe,i + /3iD- ^ ' 

It follows from (??) that: 

(1) If the dilution rate is increased at fixed feed concentrations, the steady 
state enzyme level passes through a maximum. However, at every dilution 
rate, the mixed-substrate enzyme level is always lower than the single- 
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Figure 14. Variation of the peripheral enzyme levels: (a) Activity of /3-galactosidase 
during growth of E. coli K12 on mixtures of lactose (ImM) phis ghicose (2mM), 
glucose-6-phosphate (2mM), or glycerol (4mM) (?, Fig. 3). (b) Activity of gen- 
tisate-l,2-dioxygenase and glucose-6-phosphate dehydrogenase (G6PDH) during 
growth of P. aeruginosa on saHcylate + glucose at D = 0.02 h"^ (?, Fig. 3). 
(c) Activities of alcohol oxidase and formaldehyde dehydrogenase during growth 
of H. polymorpha on various mixtures of glucose and methanol at D = 0.15 h~^ and 
Sf^t = 5 g (?, Figs. 7A,C). (d) Activity of /3-galactosidase during growth of E. 
coli K12 on various mixtures of lactose and glucose at D = 0.4 h~^ (calculated from 
Fig 5 of ?). 

substrate enzyme level. 
(2) If the fraction of Si in the feed is increased at a fixed dilution rate, the 
activity of Ei increases with The increase is linear if D is small, and 
hyperbolic if D is large. 

Both results are consistent with the data (Fig. ??). 

Based on studies wdth batch cultures of E. coli, it is wddely believed that 
glucose and glucose-6-phosphate (G6P) are strong inhibitors of /3-galactosidase 
synthesis, whereas glycerol is a weak inhibitor. It is therefore remarkable that 
the /3-galactosidase activity at any given dilution rate is the same in three 
different continuous cultures fed wdth lactose + glucose, lactose ■ G6P, and 
lactose I glycerol (Fig. ??a). The model provides a simple explanation for 
this data. To see this, let 5*1 and S2 denote lactose and glucose/G6P/glycerol, 
respectively. Now, (??) implies that any given dilution rate, the influence of 5*2 
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Figure 15. Upper panel: Variation of the glucose and galactose concentrations with 
the dilution rate during single- and mixed-substrate growth of E. coli ML308 on glu- 
cose or/and galactose (??, Appendix B). The residual concentrations of (a) glucose 
and (b) galactose during mixed-substrate growth (a) are lower than the correspond- 
ing levels during single-substrate growth (A). Lower panel: The relationship be- 
tween the substrate concentration and the peripheral enzyme level, (a) Growth of C. 
boidinii on various mixtures of methanol and glucose at D = 0.14 h~^ and Sf^t = 5.0 
g L^^ (?, Fig. 7). When the activity of alcohol oxidase increases (0 < Pmet ^ 0.4), 
the methanol concentration is constant. When the activity of alcohol oxidase satu- 
rates (/3met ^ 0.4), the methanol concentration increases, (b) Residual concentra- 
tions of galactose (open symbols) and glucose (closed symbols) during growth on 
E. coli ML308 on various mixtures of galactose and glucose at D = 0.3 h^^ and 
Sf,t = 1 (A, A), 10 (0,4), and 100 (□) mg L'^ (??). 



on the /3-galactosidase activity, Ci, is completely determined by the parameter. 

Pi ~ T7 —TP , 

YlSf,i + Y2Sf,2 

which is identical to ipi, the mass fraction of lactose in the feed, because the 
yields on lactose, glucose, G6P, and glycerol are similar. It turns out that 
the feed concentrations used in the experiments were such that ibi ~ 0.45 in 
all the experiments. Consequently, the /3-galactosidase activity at any D is 
independent of the chemical identity of 82- 
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The steady state substrate concentrations are given by the expression 



which follows immediately from (??) and (??)• Evidently, if the dilution rate is 
increased at fixed feed concentrations, the substrate concentrations increase. 
However, comparison with (??) shows that at every dilution rate, the sub- 
strate concentrations during mixed-substrate growth are lower than the cor- 
responding levels during single-substrate growth. This agrees with the data 
(Figs. ??a,b). Eq. (??) also implies that if the fraction of Si in the feed, and 
hence, (3i, is increased at a fixed dilution rate, the substrate concentration 
is constant when /3j is small, and increases when /3j is large. This occurs be- 
cause at small (3i, the linear increase of rs,i is driven entirely by the linear 
increase of Cj. Since the enzyme level saturates at sufficiently large Pi, further 
improvement of j is obtained by a corresponding increase of the substrate 
concentration. Fig ??c shows that this conclusion is consistent with the data. 
To be sure, there are instances in which the substrate concentrations appear 
to increase for all < /3i < 1 (Fig. ??d). The model implies that this occurs 
because YJK ,;A'e,i is so small that the peripheral enzymes are saturated at 
relatively small values of PiD. 

In general, one expects the mixed-substrate steady states to depend on three 
independent parameters, namely, D, (jji, and (Ty2- However, the model im- 
plies that at sufficiently small dilution rates (such that Sj ^ s/,j), the steady 
state specific uptake rates, and hence, the steady state enzyme levels and 
substrate concentrations, are completely determined by only two indepen- 
dent parameters, D and l3i{ipi). Experiments provide direct evidence support- 
ing this conclusion. Indeed, the data in Fig. ??d were obtained by perform- 
ing experiments with three different values of the total feed concentration 
{sf^t = 1, 10, 100 mg/L). It was observed that at any given fraction of galac- 
tose in the feed, the residual sugar concentrations were the same, regardless of 
the total feed concentration. Thus, the residual substrate concentrations are 
completely determined by the dilution rate and the fraction of the substrates 
in the feed (rather than the absolute values of the feed concentrations) . The 
same is true of the enzyme levels (Fig. ??b). Athough three different substrates 
(glucose, G6P, glycerol) were used in the experiments, the /5-galactosidase ac- 
tivity at any given D was the same regardless of the identity of the substrate, 
since the fraction of lactose in the feed was identical (~0.45) in all three ex- 
periments. 



3.3.2.4 High dilution rates Eq. (??) implies that if the dilution is 
increased at fixed feed concentrations, the substrate concentrations increase 
monotonically. At a sufficiently high dilution rate, S2 approaches saturating 




(54) 



40 



